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Theories that attempt to explain the observed cosmic acceleration by modifying general relativity 
all introduce a new scalar degree of freedom that is active on large scales, but is screened on small 
scales to match experiments. We demonstrate that if such screening occurs via the chameleon 
mechanism, such as in /(-/?) theory, it is possible to have order unity violation of the equivalence 
principle, despite the absence of explicit violation in the microscopic action. Namely, extended 
objects such as galaxies or constituents thereof do not all fall at the same rate. The chameleon 
mechanism can screen the scalar charge for large objects but not for small ones (large/small is 
defined by the depth of the gravitational potential, and is controlled by the scalar coupling). This 
leads to order one fluctuations in the ratio of the inertial mass to gravitational mass. We provide 
derivations in both Einstein and Jordan frames. In Jordan frame, it is no longer true that all objects 
move on geodesies; only unscreened ones, such as test particles, do. In contrast, if the scalar screening 
occurs via strong coupling, such as in the DGP braneworld model, equivalence principle violation 
occurs at a much reduced level. We propose several observational tests of the chameleon mechanism: 
1. small galaxies should accelerate faster than large galaxies, even in environments where dynamical 
friction is negligible; 2. voids defined by small galaxies would appear larger compared to standard 
expectations; 3. stars and diffuse gas in small galaxies should have different velocities, even if they 
are on the same orbits; 4. lensing and dynamical mass estimates should agree for large galaxies but 
disagree for small ones. We discuss possible pitfalls in some of these tests. The cleanest is the third 
one where mass estimate from HI rotational velocity could exceed that from stars by 30% or more. 
To avoid blanket screening of all objects, the most promising place to look is in voids. 
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I. INTRODUCTION 

The surprising finding of cosmic acceleration about ten 
years ago has motivated a number of attempts to mod- 
ify general relativity (GR) on large scales. They fall 
roughly into two classes. One involves adding curva- 
ture invariants to the Einstein-Hilbert action, the most 
popular example of which is f{R) theory [l|, [3, H, 
The other involves giving the graviton a mass or a res- 
onance width, the most well known example of which is 
the Dvali-Gabadadze-Porrati (DGP) braneworld model 
[^,@]. The first class of theories is equivalent to the classic 
scalar-tensor theories Q . The second class of theories 
is more subtle but also reduces to scalar-tensor theories 
on sub-Hubble scales In the context of DGP, 

the extra scalar can be thought of as a brane-bending 
mode. In fact, it appears the extra dimension is not even 
necessary — four dimensional ghost-free generalizations of 
the DGP model exist where it is the nonlinear interaction 
of the scalar that causes self-acceleration (TT| . 

All examples to date, therefore, introduce a light scalar 
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(or even multiple light fields) to modify gravity ^ . While 
such a scalar is welcome phenomenologically on cosmo- 
logical scales, it must be suppressed or screened on small 
scales to satisfy stringent constraints from solar system 
and terrestrial experiments. The two classes of theories 
have quite distinctive screening mechanisms. Curvature- 
invariant theories such as f{R) screen the scalar via 
the so-called chameleon mechanism, essentially by giv- 
ing the scalar a mass that depends on the local density 
[13, [H, S mi US, m HI . Theories such as DGP screen 
the scalar a la Vainshtein [2^, namely strong coupling 
effects suppress the scalar on small scales [1, . 

Our goal in this paper is to study how extended objects 
move in the presence of these two screening mechanisms. 
We will keep our discussion as general as possible even 



This is not surprising given a famous theorem due to Weinberg 
which states that a Lorentz invariant theory of a massless 
spin two field must be equivalent to GR in the low energy limit. 
Low energy modifications of gravity therefore necessarily involve 
introducing new degrees of freedom, such as a scalar. Giving the 
graviton a small mass or resonance width turns out to introduce 
a scalar d la Stueckelberg [T3l |. Theories that invoke a Lorentz 
violating massive graviton also involve a scalar, namely a ghost 
condensate [l4l llSj . Recent proposals for degravitation are yet 
another example of scalar-tensor theories in the appropriate lim- 
its [lEIil- 
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as we use f{R) and DGP as illustrative examples. In- 
deed, neither theory is completely satisfactory: accept- 
able /(i?) models accelerate the universe mostly by a 
cosmological constant rather than by a genuine modi- 
fied gravity effect [26] ; also from an effective field theory 
standpoint they are not favored in any sense over more 
general scalar-tensor theories; DGP and generalizations 
thereof [111 probably lack a relativistic ultra-violet com- 
However, it is likely the screening mecha- 
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nism, chameleon or strong coupling, will remain relevant 
in future improvements of these theories. 

The problem of motion of extended objects has a vener- 
able history. Newton famously showed in Principia that 
the Earth's motion around the Sun can be computed by 
considering the Earth as a point mass. It is worth em- 
phasizing that such a result is by no means guaranteed 
even in Newtonian gravity. This result holds only if there 
is a separation of scales, that the dimension of the Earth 
is small compared with the scale on which the Sun's grav- 
itational field varies, in other words that tides from the 
Sun can be ignored. More concretely, the Earth effec- 
tively sees a linear gradient field from the Sun; second 
derivatives of the Sun's gravitational potential can be 
safely neglected when computing the Earth's motion. 

We will work within the same zero-tide approximation 
in this paper. It has long been known that within GR, 
under the same approximation, extended objects move 
just like infinitesimal test particles. In other words, ne- 
glecting tidal effects all objects move on geodesies, or 
equivalently, their inertial mass and gravitational mass 
are exactly equal. This has been referred to as the efface- 
ment principle, that the internal structure of an object, 
no matter how complicated, has no bearing on its overall 
motion [28l |. 

It is also well known that the same is not true in scalar- 
tensor theories. In such theories, despite the universal 
coupling of the scalar to all elementary matter fields, ob- 
jects with different amounts of gravitational binding en- 
ergy will move differently. Effectively, they have different 
ratios of the inertial mass to gravitational mass, known 
as the Nordtvedt effect [29l|.^ This apparent violation 
of the equivalence principle ^ is generally small because 
most objects have only a small fraction of their total mass 
coming from their gravitational binding energy. In other 
words, such violation is suppressed by <I> where <I> is the 



^ By comparing the free fall of the Earth and the moon towards the 
sun, lunar laser ranging data already place limits on such strong 
equivalence principle violation at the level of 10~* |30| . with a 
factor of ten improvement expected in the near future [3l| . In 
this paper, we are interested in theories that respect these solar 
system bounds while exhibiting 0(1) deviation from GR in other 
environments. 

^ In this paper, equivalence principle violation means no more and 
no less than this: that different objects can fall at different rates, 
or equivalently, t hey have different ratios of inertial to gravita- 
tional mass. See [32l | for a discussion of different formulations of 
the equivalence principle. 



gravitational potential depth of the object. In the par- 
lance of post-Newtonian expansion, such violation is an 
order 1/c^ effect (c is the speed of light). Order one viola- 
tion of the equivalence principle is however observable in 
the most strongly bound objects i.e. black holes, which 
do not couple to the scalar at all — they have no scalar 
hair. This means black holes and stars fall atappreciably 
different rates in scalar-tensor theories fss!, 'sif . 

In this paper, we wish to point out that there could be 
order one violation of the equivalence principle, even for 
objects that are not strongly bound, such as galaxies and 
most of their constituents (|$| <C 1). In the parlance of 
post-Newtonian expansion, this is an order 1 effect (i.e. 
this is not even posi-Newtonian!). That this is possible 
is thanks to the new twist introduced by recent variants 
of scalar-tensor theories, namely their screening mecha- 
nisms. We will show that the chameleon mechanism leads 
to order one violation of the equivalence principle even 
for small |$| objects, while the Vainshtein mechanism 
does not. These screening mechanisms were absent in 
traditional treatments of scalar-tensor theories because 
the latter typically did not have a non-trivial scalar po- 
tential (needed for chameleon) or higher derivative in- 
teractions (needed for Vainshtein) i.e. they were of the 
Jordan-Fierz-Brans-Dicke variety [s^ [s^ |33] . 

How such equivalence principle violation comes about 
is easiest to see in Einstein frame, where the scalar me- 
diates an actual fifth force. Under the chameleon mech- 
anism, screened objects effectively have a much smaller 
scalar charge compared to unscreened objects of the same 
mass. Hence, they respond differently to the fifth force 
from other objects in the same environment; they then 
move differently. That such a gross violation of the equiv- 
alence principle is possible might come as a surprise in 
Jordan frame. One is accustomed to thinking that all ob- 
jects move on geodesies in Jordan frame. One of our main 
goals is to show that this is in fact not true. Whether 
two different objects fall at the same rate under gravity 
is of course a frame-independent issue; nonetheless, some 
readers might find a Jordan frame derivation illuminat- 
ing. Here, we follow the elegant methodology developed 
by Einstein, Infeld & Hoffmann [s^ and Damour |28j 
and apply it to situations where the screening mecha- 
nisms operate. What we hope to accomplish in part is 
putting the arguments of Khoury and Weltman [3, 1131 
on a firmer footing, who first pointed out the existence of 
the chameleon mechanism, and used Einstein frame argu- 
ments to propose equivalence principle tests using artifi- 
cial satellites. We also wish to establish that if screening 
occurs via the Vainshtein mechanism, there is no analo- 
gous order one violation of equivalence principle. 

It is worth emphasizing the different physical origins 
for equivalence principle violations of the post-Newtonian 
0(l/c^) type and the Newtonian 0(1) type. In any 
scalar-tensor theory, the scalar is always sourced by 
the trace of the matter energy-momentum tensor. In 
other words, at the level of the Einstein frame action, 
the coupling between the scalar ip and matter is of the 
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form (pT"^, where is the trace of the matter energy- 
momentum. Fundamentally, equivalence principle viola- 
tions occur whenever different objects couple to dif- 
ferently. In a theory like Brans-Dicke, this happens be- 
cause different objects of the same mass have different 
T™. The degree to which their T™ differs is of order 
(i.e. ^ the gravitational potential), and so equivalence 
principle violation of the Brans-Dicke type is typically 
small, except for relativistic objects whose T™ is highly 
suppressed i.e. most of their mass comes from gravita- 
tional binding energy. An extreme example is the black 
hole whose mass is completely gravitational in origin. In 
theories where the chameleon mechanism operates, there 
is an additional and larger source of equivalence princi- 
ple violation, and that is through variations in the effec- 
tive scalar charge that a macroscopic object carries. The 
scalar self-interaction in a chameleon theory makes the 
scalar charge sensitive to the internal structure of the 
object; the scalar charge is not simply given by a vol- 
ume integral of T™ as in Brans-Dicke, but is given by a 
screened version of the same. As we will see, chameleon 
screening could introduce 0(1) fluctuations in the scalar 
charge to mass ratio, even for non-relativistic objects. 
What is perhaps surprising is that this happens despite 
an explicitly universal scalar coupling at the level of the 
microscopic action. One can view this as a classical renor- 
malization of the scalar charge. We will also see that in 
theories where the Vainshtein mechanism operates, there 
is no such 0(1) charge renormalization. Therefore these 
theories only have equivalence principle violation of the 
Brans-Dicke, or post-Newtonian, type. 

Equivalence principle violations of both types dis- 
cussed above, Newtonian or post-Newtonian, should be 
clearly distinguished from yet another type of violation, 
which is most commonly discussed in the cosmology lit- 
erature — where there is explicit breaking of equivalence 
principle at the level of the microscopic action, e.g. ele- 
mentary baryons and dark matter particles are coupled to 
the scalar with different strengths [H, HO, [U, IH, [ll H] . 
In this paper, we are interested in theories where this 
breaking is not put in 'by hand' i.e. it does not exist 
at the level of the microscopic action; rather, the equiv- 
alence principle violation comes about for macroscopic 
objects by way of differences in response to the screening 
mechanism. 

Further experimental tests come to mind after our 
investigation of the problem of motion under the two 
screening mechanisms. We propose several observational 
tests using either the bulk motion of galaxies or the in- 
ternal motion of objects within galaxies. The chameleon 
mechanism will predict gross violations of the equivalence 
principle for suitably chosen galaxies while the Vainshtein 
mechanism will not. Care must be taken in carrying out 
some of these tests, however, to avoid possible Yukawa 
suppression and/or confusion with astrophysical effects. 
We point out that the most robust test probably comes 
from the measurement of internal stellar and HI velocities 
of small galaxies in voids. Systematic differences between 



the two, or lack thereof, will put interesting constraints 
or rule out some of the current modified gravity theories. 

A guide for the readers. Readers who are interested 
primarily in observational tests can skip directly to WII 
where relevant results from the more theoretical earlier 
sections are summarized. ^TUserves as a warm-up exercise 
wherein we derive the equality of inertial mass and grav- 
itational mass in GR for extended objects. ^IIII demvsti- 
fies this surprising effacement property of GR by viewing 
GR through the prism of a field theory of massless spin 
two particles. This section is somewhat of a digression 
and can be skipped by readers not interested in the field 
theoretic perspective. ijlVl is where we deduce the motion 
of extended objects under the chameleon mechanism in 
Einstein frame. The Jordan frame derivation is given in 
the Appendix. In fjV]we do the same for the Vainshtein 
or strong-coupling mechanism. Observational tests are 
discussed in WII where we also summarize our results. 

A number of technical details are relegated to the Ap- 
pendices. Appendix discusses how the requisite surface 
integral in the Einstein-Infeld-Hoffmann approach can be 
done very close to the object. Appendix [B] maps our 
parameterization of general scalar-tensor theories to the 
special cases of Brans-Dicke and f{R). Appendix [C] pro- 
vides the Jordan frame derivation. Appendix [D] derives 
in detail the equation of motion under the Vainshtein 
mechanism. Appendix |E] supplies useful results for an 
expanding universe. 

A brief word on our notation: the speed of light c is 
set to unity unless otherwise stated; 8ttG and 1/Afpj are 
used interchangeably; the scalar field ip (or tt) is chosen 
to be dimensionless rather than have the dimension of 
mass. 



II. THE PROBLEM OF MOTION IN GENERAL 
RELATIVITY — A WARM UP EXERCISE 

We start with a review of an approach to the prob- 
lem of motion introduced by Einstein, Infeld & Hoffmann 
[sst and developed by Damour [11] , using GR as a warm 
up exercise. The idea is to work out the motion of an 
extended object using energy momentum conservation, 
rather than to assume geodesic motion from the begin- 
ning. It helps clarify the notion that the motion of ex- 
tended objects is something that should be derived rather 
than assumed — depending on one's theory of gravity, the 
extended objects might or might not move like infinites- 
imal test particles. 

The Einstein equation G^'' = SttG T'"^'' can be rewrit- 
ten as: 

GW/ = 87rGt/ (1) 

where G'-^-'^'' is the part of the Einstein tensor G^'^ that is 
first order in metric perturbations, and i^'' is the pseudo 
energy-momentum tensor which is related to the energy- 
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momentum tensor T™^^ 



by 



(2) 



where G^'^\" is the higher order part of G^^ , including 
but not necessarily limited to second order terms. The 
superscript m on T™^" is there to remind us this is the 
matter energy-momentum, in anticipation of other kinds 
of energy-momentum we will encounter later. Note that 
the identity d^O'^'^^^ = implies t^" is conserved in the 
ordinary (as opposed to covariant) sense. 

Let us draw a sphere of radius r enclosing the extended 
object whose motion we are interested in. The linear 
momentum in the «-th direction is given by 



F, = / (Px t 



(3) 



where the volume integral is over the interior of the 
sphere. We assume that ti' is dominated by the extended 
object such that Pi is an accurate measure of the linear 
momentum of the object itself, even though we are inte- 
grating over a sphere that is larger than the object. 
The gravitational force on the extended object is 



d^xdoU" 



d^xdjU^ 



(4) 



where in the last equality we have used the Gauss law 
to convert a volume integral to a surface integral i.e. 
dSj — dAxj, where dA is a surface area element and 
X is the unit outward normal. As expected, the gravi- 
tational force is given by the integrated momentum flux 
through the surface. This is the essence of the method 
introduced by Einstein, Infeld & Hoffmann: the problem 
of motion is reduced to performing a surface integral — 
i.e. there is no need to worry about internal forces within 
the object of interest. 

Assuming that T"^i^ at the surface of the sphere is 
negligible, we need only consider the contribution to 
from G(2),J' (Eq. @). And if the surface is located at 
a sufficiently large distance away from the object of in- 
terest, we can assume the metric perturbations are small 
and only second order contributions to G^^-'i-' need to be 
considered (without assuming metric perturbations are 
small close to the object). In the Newtonian gauge. 



ds^ = -(1 + 2^)df + (1 - 2^^)6^Jdx'dx^ 



(5) 



we have: 



-di^dj'^ - ^^'^^J<^> + 2^G'^^\^ (6) 



ignoring time derivatives. On the right hand side, we 
are being cavalier about the placement of upper/lower 
indices — it does not matter as long as we are keeping 



only second order terms. The symbol denotes dkd'^. 
The first order Einstein tensor G'^^\^ is given by: 



(7) 



again, ignoring time derivatives. 

The next step is to decompose $ and ^ around the 
surface of the sphere into two parts: 



$ = $0 + '^i{r) 



(8) 



where $o and represent the large scale fields due to 
the environment/background while $1 and ^1 are the 
fields due to the extended object itself. This decomposi- 
tion can be formally defined as follows: $1 and solve 
the Einstein equation with the extended object as the 
one and only source; $0 and are linear gradient fields 
that can always be added to solutions of such an equation 
(because the equation consists of second derivatives) . We 
have in mind $0 and 'So are generated by other sources 
in the environment, and they vary gently on the scale of 
the spherical surface that encloses the object i.e. their 
second gradients can be ignored: 



$0(2^) ^ $o(0)-Hft$o(0) 
*o(2?) -*o(0)+ 9,^-0(0) 



(9) 



where denotes the center of the sphere. In other words, 
we assume di^Q and 9i^o hardly vary on the scale of 
the sphere. On the other hand $1 and 5*1, the object 
generated fields, have large variations within the sphere. 
At a sufficiently large r, we expect both to be 



-GM/r 



(10) 



i.e. dominated by the monopole. 

Note that we have carefully chosen the radius r of the 
enclosing sphere to fulfill two different requirements: that 
r is smaller than the scale of variation of the background 
fields, and that r is sufficiently large to have the monopole 
dominate. The latter condition simplifies our calculation 
but is not necessary in fact, as is shown in Appendix [XI 
Also we have in mind a situation where the density within 
the object of interest is much higher than its immediate 
environment — that is the definition of 'object' after all — 
so that the total mass inside our sphere is dominated by 
the object's. 

Plugging the decomposition (Eq. ([S])) into Eq. © and 
performing the surface integral as in Eq. ([5]), we find ^ 



1 



SttG 



dSjG^^\^ 
3^ 



3 dr 



(11) 



* We have not assumed equality of <I>o ^^nd 'I'Oi nor <I>i and 'I'l in 
Eq. mil l. Keeping their relationship general will be useful for 
the scalar-tensor case. Eq. JTTJ is derived assuming only that 
# and ^ are subject to an object-background split with a linear 
background. 



5 



This result has several remarkable features. First, only 
cross-terms between the background fields (subscript 0) 
and the object fields (subscript 1) survive the surface 
integration. A term such as ^ dSjdi^odj'^'o vanishes 
because by assumption, (?i$o and dj'^o are both con- 
stant (or nearly so on the scale of the sphere). A term 
such as f dSjdi^idj'i/i vanishes because both $i and 
^'i are spherically symmetric. Second, only terms pro- 
portional to di^o remain. There are no di'^o terms. This 
should come as no surprise — by ignoring time deriva- 
tives, we are in essence assuming the object is mov- 
ing with non-relativistic speed and so its motion should 
only be sensitive to the time-time part of the back- 
ground/environmental metric. 

Finally, plugging the monopole <i>i and 5*1 into Eq. 
(fTTj) gives us: 



(12) 



which is the expected GR prediction in the Newtonian 
limit. 

To complete this discussion, it is useful to note that 
the mass of the object is given by 



M = 



(13) 



Its time derivative AI = J d^x dito" can once again be 
turned into a surface integral M — § dSt t^" . Assuming 
the energy flux through the surface is small, we can ap- 
proximate M as constant. The center of mass coordinate 
of the object is defined by 



X X 



(14) 



Making the constant M approximation, its time deriva- 
tive is given by 



d^xdj{x'to^)-W 



/M 



d^xto'/M . 



(15) 

This is precisely Pi /M defined in Eq. ^ . Taking another 
time derivative, we can see that Eq. (fT2)) is equivalent to 



MX' 



(16) 



As expected, the center of mass of an extended object 
moves on a geodesic. Its mass cancels out of the equa- 
tion of motion i.e. inertial mass equals gravitational mass. 
That this holds independent of the internal structure of 
the extended object is sometimes referred to as the ef- 
facement principle ^2811. 

It is useful to summarize the approximations or as- 
sumptions we have made to obtain this result, some of 
which are crucial, some of which not. 

1. We have taken the non-relativistic/Newtonian limit 
i.e. the object has the equation of state of dust and is 
slowly moving, and time derivatives of the metric are 



ignored. These assumptions can be relaxed and the ob- 
ject would still move on geodesies. For instance, photons 
surely move on geodesies in GR, but we will not dwell on 
the proof here. 

2. We assume the metric perturbations are small at 
the surface of the sphere that encloses the object of inter- 
est. The derivation does not assume however that metric 
perturbations are small close to the object. For instance, 
Eq. PB)) remains valid for the motion of a black hole. 
Note also that the assumption of a spherical surface is 
not crucial — the surface can take any shape, as discussed 
in Appendix [Xl 

3. We assume there is a separation of scales, that ex- 
ternal tides can be ignored on the scale r of the enclosing 
surface, and that the monopole contribution dominates 
at r. The last assumption can be easily relaxed. As 
shown in Appendix]^ one could choose an r that is close 
to the object, and the higher multipoles would still not 
contribute to the relevant surface integral. This is an im- 
portant point because the enclosing surface is purely a 
mathematical device — it should not matter what shape 
or size it takes, as long as it is not so big that too much 
of the exterior energy-momentum is enclosed or, equiva- 
lently, external tides become non-negligible. 

Let us end this section with an interesting observa- 
tion. This derivation makes use of the metric that the 
extended object sources (Eq. pH)) ). It is intriguing that 
what the object sources also determines how the object 
responds to an external background field. In fact this is 
not surprising from the Lagrangian viewpoint, and is in 
a sense the modern generalization of Newton's third law: 
the same object-field interaction Lagrangian term yields 
both the source term in the field equations for gravity and 
the gravitational force term in the object's equation of 
motion. This will be useful in understanding how equiv- 
alence principle violation comes about in scalar-tensor 
theories. 



III. FIELD-THEORETICAL CONSIDERATIONS 

The fact that when tidal effects are negligible, ex- 
tended objects move along geodesies of the background 
metric, looks somewhat magical if derived as above but 
is in fact forced upon us by Lorentz invariance.^ In par- 
ticular, it extends beyond the non-relativistic/ weak- field 
regime discussed above. For example, a small black-hole 
falling in a very long-wavelength external gravitational 
field, will follow a geodesic of the latter. 

This can be shown by imposing the requirement that 
the theory of massless spin-two particles (gravitons) in- 



In this section, we digress a bit and wish to throw some light 
on the effacement property of GR from a field-theoretical per- 
spective. Readers interested mainly in the problem of motion in 
scalar-tensor theories can skip to the next section without loss 
in the continuity of the logic. 
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teracting with matter be Lorentz invariant. In this sec- 
tion, we will thus abandon momentarily the geometric de- 
scription of GR, and treat gravity as a theory of gravitons 
in Minkowski space. After all, from the modern quantum 
field theory (QFT) viewpoint this approach is what legit- 
imizes GR as a theory of gravity in the first place: GR is 
the only consistent low-energy Lorentz-invariant theory 
of interacting massless spin-two particles. 

So, let's describe gravitons in flat space via a tensor 
field /i^j^. We want to exploit the smallness of the ob- 
ject in question with respect to the typical variation- 
scale of the external gravitational field. What we can 
do then, is treat the object as a point-particle with tra- 
jectory x'^(t), and couple it to the graviton field in all 
possible ways allowed by symmetry. (Here the x^^s are 
flat-space coordinates, and t is flat-space proper-time: 
dr^ — —rj^iy dx^dx'^ .) This yields an effective field the- 
ory valid for wavelengths much longer than the object's 
size. The important symmetry for us is just Lorentz- 
invariance. When applied to gravitons, this also demands 
invariance under gauge transformations 

hf^„ hf_,^ + a^^p + df_,^t, , (17) 

see e.g. ref. j45| . However at the order we will be working, 
this is not relevant. The low-energy effective action is 

S'cff = ~M J dr + SgrsiA^T't^i^] 

-Mf JdTh^.S'^'' + ... (18) 

The first line is just the sum of the free point-particle 
action and of the pure gravity action, which includes 
the graviton free action as well as all graviton self- 
interactions. In the second line, we explicitly display one 
interaction term between h^^^, and the point-particle. S'^'^ 
is a Lorentz-tensor built out of the point-particle degrees 
of freedom: its trajectory x^{t), as well as possibly inter- 
nal degrees of freedom such as spin, multipole moments, 
etc.; / is a coupling constant. For convenience, we pulled 
out of / an explicit factor of the particle's mass M. The 
dots stand for additional interactions, which involve ei- 
ther higher powers of /i^j/, or derivatives thereof, or both. 

Now, Weinberg (46l| showed that for this effective the- 
ory to be Lorentz-invariant, in the low-energy graviton 
limit: 

i) S^'^ must be totally independent of the internal de- 
grees of freedom, and given by 

5'"^ = i^i"^ . (19) 



redefinition of h^^. Let's therefore set it to 1/2, so that 
the effective action becomes 

S-eff = ~M JdT{l + ^h^.X^x") + Sgr..Ku] + ■■■ (20) 

We recognize in parentheses the first order expansion of 
the point-particle action in GR: 



^ „A/ J dT^-g^,x>^x- (21) 

with metric g^i^ = rj^^ -f h^^. In fact, one can show 
that because of the gauge-invariance (fT7| (which is also a 
consequence of Lorentz-invariance), the non-linear inter- 
actions between hf^i, and the point-particle (which make 
up part of the 'dots' in eq. ([20]) ). as well as the (low- 
energy) graviton self- interaction of S'grav, must combine 
to yield the usual GR expressions. The remainder of the 
'dots' — interactions involving higher derivatives of h^^, — 
encode non-minimal couplings of h^^ with the particle's 
multipole moments f47j, which lead to tidal effects. These 
however are negligible for very long wavelength gravita- 
tional fields, their relative importance being suppressed 
by powers of a/A, where a is the object's size, and A the 
graviton's wavelength. 

These are of course very well known notions. We re- 
viewed them here, to contrast the very constrained situ- 
ation of GR with a generic scalar-tensor theory. For GR, 
any object will source and respond to long-wavelength 
gravitational fields as if it were a point particle: the 
equivalence between inertial and gravitational masses is 
guaranteed by Lorentz invariance, and as a consequence 
it is totally independent of the object's internal struc- 
ture, and of the forces that are keeping the object intact 
(including the gravitational one). It applies equally to 
macroscopic objects with non-relativistic structure (like 
e.g. a very soft, very cold gas cloud) and to those with 
a very relativistic one (like e.g. a black hole). Notice 
that the parameter M appearing in eq. (pO| is the total 
physical mass of the object, which includes all classical 
and quantum contributions to it. Among the former we 
have, on top of the constituents' masses, the classical 
binding energy of a macroscopic system. In fact this is 
just the tree-level approximation of the quantum contri- 
butions, which are particularly relevant for microscopic 
systems where the bulk of the physical mass is quantum- 
mechanical in origin, like e.g. for the proton. Wein- 
berg's argument [4^ shows that upon quantum correc- 
tions, the inertial and gravitational masses of a quantum- 
mechanical system get renormalized by exactly the same 
multiplicative factor ^. 



ii) The coupling constant / must be universal, that is, 
the same for all particles. 

This is nothing but the equivalence "principle" , which 
from this viewpoint is rather a theorem. Indeed, since 
/ is a universal constant, we can absorb its value into a 



This is tlie gravitational analogue of a similar — but not quite as 
powerful — statement for electric charges: quantum corrections 
renormalize all charges by the same multiplicative factor, be- 
cause charge renormalization can be thought of as wave-function 
renormalization for the photon field. See e.g. ref. |45ll . 
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No analogous statement can be made for scalar cou- 
plings. Consider indeed a scalar-tensor theory for gravity. 
Also assume that it is a 'metric' theory, which means that 
there is a conformal frame (Jordan's) where matter-fields 
are all minimally coupled to the same metric. This is the 
closest we can get to implementing the equivalence prin- 
ciple. Now, upon demixing the scalar Lp from the graviton 
h^u (that is, in Einstein frame), this prescription corre- 
sponds to a direct coupling </? T™ between the scalar and 
matter fields, where = T™^'^ is the trace of the mat- 
ter energy-momentum tensor. The problem with such a 
coupling is that T™ is not associated with any conserved 
charges. That is, in the above point-particle approxima- 
tion, extending GR to a scalar-tensor theory (in Einstein 
frame) amounts to adding an interaction term 



Sint = Q I drip 



(22) 



as well as of course a Lagrangian for Lp, its self- 
interactions, and its interactions with /i^^. The dots in 
Siat denote interactions involving derivatives of ip, which 
can be thought of as multipole interactions — negligible 
for a long-wavelength scalar field. Now, the coupling Q 
is the spatial integral of T™ over the object's volume, 
and is not a conserved charge. So, for instance, if we 
consider a collapsing sphere of dust eventually forming 
a black-hole, its coupling to h^^ (the mass M) is con- 
served throughout the collapse, whereas its coupling to 
If (the scalar charge Q) is not: as the dust particles be- 
come more and more relativistic, their T™ approaches 
zero, even as the total mass M remains constant, domi- 
nated more and more by gravitational binding energy. In 
fact we already know that for a black-hole, which is the 
end product of our example, Q must be zero, because of 
the no-hair theorem [i^: a black-hole cannot sustain a 
scalar profile, which in our language means that in the 
point-particle approximation it cannot couple to a long- 
wavelength ip. So, in comparing how a black-hole and a 
non-relativistic gas cloud fall in a very long-wavelength 
scalar-tensor gravitational field, we expect order one vi- 
olation of the equivalence principle: the black-hole will 
not follow a geodesic of the Jordan-frame metric. 

In general, even if we start with a universal coupling 
(fT"^ at the microscopic level, the ratio Q/M for ex- 
tended objects will be object dependent, and sensitive 
to how relativistic the object's internal structure is: it 
will be one for very non-relativistic objects (with a suit- 
able normalization of i^), and zero for black holes. This 
kind of variation in Q/M , which is relativistic in nature 
(0(l/c^) in a post-Newtonian expansion), is available to 
all scalar-tensor theories. As we will see, theories where 
the chameleon mechanism is at work possess an addi- 
tional and more dramatic source of Q/M variation — it 
can be 0(1) even for non-relativistic objects. This screen- 
ing mechanism arises from a particular kind of scalar self- 
interaction, and it results in a scalar charge Q not just 
determined by an integral of T™, but also by the precise 
non-linear dynamics of the scalar itself close to the ob- 



ject. Essentially, the scalar field profile sourced by the 
object at short distances contributes, via non-linearities, 
to the coupling of the object to a long-wavelength exter- 
nal scalar. (As we will see, for the Vainshtein mechanism 
this scalar charge renormalization does not occur.) 

The lesson is this: from the low-energy point-particle 
action viewpoint, the charges Q for different objects are 
free parameters, and are not forced to be related to 
the corresponding total inertial masses M in any way. 
Therein lies the seed for equivalence principle violation. 



IV. THE PROBLEM OF MOTION IN 
SCALAR-TENSOR THEORIES — SCREENING 
BY THE CHAMELEON MECHANISM 

The fact that different objects might fall at different 
rates in scalar-tensor theories is easiest to see in Einstein 
frame — the scalar ip mediates a fifth force and the equiv- 
alence principle is violated whenever the scalar charge 
to mass ratio fluctuates between objects. On the other 
hand, that this should be true in Jordan frame might 
come as a bit of a surprise: after all, are not all objects 
supposed to move on geodesies in Jordan frame? We 
therefore present two different derivations, one in each 
frame. The Einstein frame derivation is presented in 
! jIV Al The Jordan frame derivation is briefly summa- 
rized in ijlVBI with details given in Appendix [Cl In addi- 
tion. Appendix [B] describes how to map our description 
of general scalar-tensor theories to the special cases of 
Brans-Dicke and f{R). 

The Einstein frame action is 



S = M^i / d-x 



ii?-iV^^VV-l/(^) 



(23) 



Here, "denotes quantities in Einstein frame, R is the Ein- 
stein frame Ricci scalar, V ^tpS/^tp is supposed to remind 
us that the contraction is performed with the Einstein 
frame metric (there is otherwise no difference be- 
tween V ^ip) and V^(p: V^(^ = "S/ ^ip = d^ip). The symbol 
ip denotes some scalar that would contribute to a fifth 
force. Upon a field-redefinition, its kinetic term can al- 
ways be written as in eq. (j23p . For notational conve- 



nience for what follows, we decide to factor an M|[ out 
of the whole scalar-tensor gravitational action. With this 
normalization Lp is dimensionless. We assume has a 
potential V{lp), which given our normalization has mass- 
dimension two. The symbol ipm denotes some matter 
field whose precise nature is not so important; the most 
important point is that it is minimally coupled not to g^y 
but to {Lp)g ^i, , which is precisely the Jordan frame 
metric: 

= f7-2(<^)g^^ . (24) 
Performing the conformal transformation gives us the 
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Jordan frame action: 



S = Af 2j / d 



-Vt'^{ip)V 



where 



3 /dlnn^ 



(25) 



(26) 



The matter energy-momentum in the two frames (defined 
as the derivative of the action w.r.t. the corresponding 
metric) are related by: 



(27) 



In anticipation of the fact that we will be carrying out 
perturbative computations, let us point out that fl'^{(p) 
must be close to unity for the metric perturbations to be 
small in both Einstein and Jordan frames. We can thus 
approximate it at linear order in ip, 



(28) 



where a is a constant, and \aip\ <ti 1. We are not assum- 
ing, however, that the fractional scalar field perturbation 
(i.e. Sip/if) is in any way small. In fact, nonlinear effects 
in Sip/ip are going to be crucial in enabling the chameleon 
mechanism. Under the same approximation, a useful re- 
lation is 



dip 



-2a 



(29) 



The reader might find it helpful to keep in mind a partic- 
ularly common (if not particularly natural) form of the 
conformal factor, 57^ = exp [— 2a(/3] with constant a, in 
which case Eq. ((29)) is exact and Eq. ((28)) is a Taylor 
expansion. 

As elaborated in Appendix [Bl a — l/\/6 in f{R) mod- 
els and wbd = (1 — 6a^)/(4a^) in Brans-Dicke theory. 



A. Derivation in Einstein Frame 



The Einstein frame metric equation is 



GfJ^ = 87rG 



(30) 



where 



^ 8^G 



^ [v^^ifVip - s^'' [iv„(^v"^ + y] } .(31) 

The scalar field equation is 



Following the derivation in the GR case f^ fll)) . we perform 
an object-background split, derive solutions to the object 
fields, and compute an integral of momentum flux over 
an enclosing sphere to find the net gravitational force on 
the extended object. 



Object-background Split 



We work in the Newtonian gauge: 



-(1 + 2^)dt^ + (1 - 2'^)S^jdx'dx' 



(33) 



As before (Eq. ([5])), we carry out a decomposition of the 
fields at some radius r away from the object: 

$ = l>o + ^i(r) , * = *o + *i(r) , ip = ipa + '/'I (0(34) 

where the subscript denotes the back- 
ground/environmental fields and the subscript 1 
denotes the object fields. The background fields are 
assumed to be well approximated by linear gradients: 



l>o(x)~l>o(0)-Ha,l>o(0)x^ 

^o(2^)-*o(o) + a,^'o(o)x' 

i^o(^) — </5* + diipQ{0)x'' 



(35) 



where denotes the origin centered at the object, and 
we use to denote the background scalar field value 
there. The assumption here is that these background 
fields sourced by the object's environment vary on a scale 
much larger than r, where we will draw our enclosing 
sphere. 

It is worth emphasizing that this object-background 
split needs only make sense at the surface of the enclos- 
ing sphere (over which we will compute the integrated 
momentum flux), but not necessarily at the location of 
the object itself. 



Solutions for tfie Object Fields 

To work out solutions for the object fields, we first 
examine the Einstein equation linearized in metric per- 
turbations: 

V^* = 47rGp - i^oV'^V + jd.ipd'ip + 



(a,a,-i^yV2) (§-$) = 



5^*-)-iv2(l> 



\{-ldopd\-\dkvd''^-W) (36) 

where we have assumed the matter is non-relativistic, 
characterized only by its energy density p. 

We can ignore all second order scalar field terms on 
the right hand side of these equations. This is because 
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as we will see shortly, Lp is already first order in G. Thus 
its contribution to Einstein's equations is on an equal 
footing with post-Newtonian corrections, which we are 
neglecting. In addition, we will ignore V as well. The 
rationale is that Gp ^ V inside the object of interest, 
whereas outside the object, V or any other sources of 
energy-momentum are just small (in the sense that the 
total mass within the enclosing sphere is dominated by 
the object). Essentially, we assume the metric in the 
vicinity of the object is sourced mainly by matter rather 
than the scalar field :^ 



Under the same assumptions as in i HIl we get 

$1 = *i = -GM/r 



(37) 



(38) 



where M is the mass of the object. 

As in the case of GR, we in essence assume there is a 
separation of scales: that r can be chosen small enough 
that second gradients of the background tidal fields can 
be ignored, but large enough for the monopole to dom- 
inate. The latter condition is helpful but not essential 
(see Appendix [X]). 

The next task is to compute the scalar field profile 
sourced by the object. The scalar field equation (Eq. 
([5^ ) can be written as 



. dV 
VV = -T^+a%T:Gp 

dip 



(39) 



where we have ignored time derivatives and corrections 
due to metric perturbations, and we have used the ap- 
proximation Eq. (|29p . Here, it is crucial we do not as- 
sume 5lp/lp is small even though </? is. The important 




Under this assumption, the anisotropic stress vanishes. This is 
emphatically an Einstein frame statement. The same will not be 
true in Jordan frame. 



FIG. 1: A scalar potential for the chameleon mechanism. 
The effective potential felt by ip in the presence of sources is 
the sum of the self-interaction potential V{ip) and the scalar- 
matter coupling (Q87rG)p</J. This can give tp a large mass at 
locations where the matter density p is large. 

nonlinearity for the chameleon mechanism is in the po- 
tential term. For instance, for the potential of fig. [T]the 
chameleon mechanism could operate, meaning that for an 
object of the right size and density, ip could be trapped at 
some small value inside much of the extended object — a 
value where the two terms on the right hand side of Eq. 
([55)1 roughly balance. This corresponds to ip having a 
large mass, or a small Compton wavelength, inside the 
object (see Fig. [1]). The object will therefore have an 
exterior scalar profile that is screened i.e. only a frac- 
tion of the object's mass (the fraction that happens to 
live within a shell approximately a Compton wavelength 
thick at the object's boundary) sources the exterior scalar 
profile. In other words, we expect the object to produce 
a scalar field that at large r is dominated by a monopole 
of the form: 



Lpi{r) = -ea 



2GM 



(40) 



where e quantifies the degree of screening. 

If the scalar field remains light inside the object, we 
have dV/dip <^ aSnGp in which case the exterior scalar 
profile is sourced by the full mass of the object, i.e. e = 1 
in Eq. (HO]). On the other hand, suppose the scalar field 
is massive and therefore short-ranged inside the object, 
the exterior scalar profile is sourced only by the mass 
residing in a shell at the boundary of the object. Sup- 
pose the object has a radius rc and the shell has thickness 
Arc- The scalar field sits at some small equilibrium value 
(where the potential and density terms on the right hand 
side of Eq. (|39p cancel) at all radii interior to the object 
- that is, except where the shell is. At the shell, the 
scalar field starts to pull away from this small expecta- 
tion value because of the force exerted by the density 
term. In other words, at the shell, Eq. ([55]) can be ap- 
proximated by 5(p/ (r^Arc) ^ a SttG M/(47rrj?/3), where 
Sip ^ ip^ — ipc — ip* (ipc is the equilibrium value where 
the the potential and density terms balance in the deep 
interior of the object, and ip^ is approximately the scalar 
field value far away from the object). We therefore have 
iAvc/rc ^ ((/3*/2a) • {GM/rc)~^. This is precisely the 
fraction by volume of the object that sources the exte- 
rior scalar field. This means the suppression factor e in 
Eq. go]) should be 



ip* 
2a 



GM 



if 



tp* 
2a 



GM 



1 -1 



< 1 



(41) 



where GM / Vc is the gravitational potential of the object. 
The inequality is required for consistency: it had better 
be true that SAr^/rc < 1 for screening to occur. This is 
known as the thin-shell condition. When this condition 
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is violated, there is no screening: 



1 if 



2^ 



GM 



> 1 



(42) 



essentially meaning that aSnGp ^ dV/dip throughout 
the object, in which case the scalar field Lp cannot even 
reach the equilibrium value (pc inside the object. It is 
interesting to note that the thin-shell condition does not 
depend explicitly on the potential V . The implicit re- 
quirement is that the potential is of such a form that 
equilibrium (between density and potential) is possible. 
A more detailed derivation of the thin-shcU condition can 
be found in [20j . 

Along the lines of ^IIW it is useful to introduce the idea 
of a scalar charge: by analogy with — —GM/r for 
the gravitational potential, let us write the scalar profile 
as ^pi = —2GQ/r (the factor of 2 is motivated by the 
appearance of 87rp rather than Anp in Eq. (IHS)) ). By 
inspecting Eq. (|40p. we can see that the dimcnsionless 
scalar charge Q is 



Q EE eaM 



(43) 



Screened objects carry a much reduced scalar charge as 
viewed from outside. 

Just as in the case of the metric fluctuations, one can 
always add to ipi a linear gradient of the form Eq. (j35p 
to represent the field generated by sources in the envi- 
ronment. To summarize, $, ^ and ip are given by Eqs. 
(131, dSSl), iHl) and (gni). 



Surface Integral of Momentum Flux 



but this can be ignored by our assumption that the scalar 
field </?i is sourced mainly by the object rather than its 
immediate environment. 

Putting together all contributions (Eq. & Einstein 
frame analog of Eq. (fTTj) ). we therefore have 



2 r 



r 
2G 



9^$o - 



4 9^1 2 9$i\ ^ dipi 
3 dr 3 dr ) dr 



(47) 



Putting our solutions for the various fields into the 
above expression, we obtain 



MX' 



= -M 



M 



(48) 



where have used Eq. and we have equated Pi and 

MX*, where is the object center of mass coordinate 
(Eq. p5|) ). This is the main result of this section. The 
fact that screened and unscreened objects have different 
scalar charge to mass ratios Q/M (or equivalently, dif- 
ferent e) means they fall at different rates. This is the 
origin of apparent equivalence principle violation. 

This result can be further simplified if we make one 
more assumption, which does not necessarily hold in gen- 
eral. Both $0 and ipo are sourced by the environment. 
Suppose they are sourced in a similar way: 



VVo 



iirG p 



(49) 



where p is the environmental density. The latter equation 
follows from Eq. (|39p provided that potential terms are 
small compared to density terms on the scales of interest. 
Eq. (|49| then implies 



We are ready to compute the gravitational force fol- 
lowing Eq. (H: 



(44) 



where the surface integral is again over some sphere en- 
closing the object of interest. By analogy to Eqs. ([1]) and 
(PI), the pseudo energy-momentum tensor is: 



I 1/ rpm V 



SttG ^ 



(45) 



As before, T™^-' outside the object is small, and we have 
already computed the contribution (just put ^ on 

top of quantities in Eq. fTT]) '). What remains to be done 
is to compute the contribution from T'^^-' (Eq. (|3ip ). 
The same tricks that work for the G^'^\^ contribution 
work here, and we obtain 



1 2 o d'Pi 
2g' 



(46) 



Note that there should be a contribution from the poten- 



tial of the form: 



G 



dSjV 



^i'inry3)d,ipodV/dip\^„ 



po 2a $0 ■ 
Putting this into Eq. (|48)l . we therefore obtain 



(50) 



MX' = -A/(9,$o 

= ^Mdi^Q [1 + 2ea^] 



Q 



(51) 



The fact that this result depends on the square of 
a makes it clear that the sign of a has no physical 
meaning — chameleon mechanism only requires its sign be 
opposite to that of dV/dip (see Eq. ([M]) ). 

Suppose one has a scalar-tensor theory where 2a^ is 
order unity, as is true in f{R) theories, then screened ob- 
jects (e ^ 0) and unscreened objects (e = 1) would move 
on very different trajectories i.e. order unity apparent vi- 
olation of the equivalence principle. Note that infinites- 
imal test particles should go unscreened and therefore 
have e = 1. 

Let us end with a summary of further approximations 
we have made in this section to augment those listed at 
the end of flU 

1'. We have taken the non-relativistic or quasi-static 
limit i.e. time derivatives of the scalar field are ignored. 
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2'. We assume the scalar field does not contribute 
much to the r.h.s. of Einstein's equations ([55]) . 

3'. We assume there is a separation of scales, that sec- 
ond derivatives of the background/environmental scalar 
field can be ignored on the scale r of the enclosing sur- 
face, and that the monopole contribution to the object's 
scalar profile dominates at r. As discussed in Appendix 
1X1 this latter assumption can be relaxed. The former 
assumption implicitly requires the Compton wavelength 
of the exterior scalar field to be large compared with r. 
The scalar field would be Yukawa suppressed if this were 
not true. This has implications for observational tests as 
we will discuss in i jVIl 

To this list, which mirrors the list in f|lll we should 
add: 

4'. We assume the conformal factor ^"^{^p) is close to 
unity, or equivalently \atf\ <C 1 (see Eq. (|28)) V This as- 
sumption is not strictly necessary in our Einstein frame 
calculation. It will become necessary when we try to 
compare the Jordan frame calculation with it. One con- 
sequence of this assumption is that the Einstein frame 
and Jordan frame masses of the object are roughly equal 



M = / (Fxp 



X p 



(52) 



n^p, recalling that T 



m 




-p and T^qO 



since p 

In summary, the main results of this section are encap- 
sulated in Eqs. (gT]), ([48]) and ([5l|), which make differ- 
ent levels of approximations, starting from the minimum. 
The key observation is that, despite the explicitly uni- 
versal coupling of the scalar to matter at the microscopic 
level (Eq. ((23| ). scalar charges of macroscopic objects 
get classically renormalized by nonlinear chameleon ef- 
fects. This produces fluctuations in the scalar charge to 
mass ratio Q/M for different objects, leading to different 
rates of gravitational free fall. Screened objects, which 
have Q/M ~ 0, accelerate less than unscreened objects. 
As we will see, no such (non-relativistic) renormalization 
occurs for the Vainshtein mechanism. 



B. Derivation in Jordan Frame 

In principle, the simplest way to obtain the Jordan 
frame results is to transform the end results of ijlV Al 
Eqs. (|17)) . and ([?T|) . into Jordan frame by a con- 
formal transformation g^^ = Vl~'^g^,y. After all, whether 
two different objects fall at different rates is a frame- 
independent issue. This might not, however, satisfy read- 
ers who are accustomed to the Jordan frame viewpoint: 
should not everything move on geodesies? In this sec- 
tion, we therefore briefly sketch a Jordan frame deriva- 
tion. Details can be found in Appendix [Cl 

The Jordan frame metric in Newtonian gauge is 



-(1 + 2^)df + (1 - 2'^)6ijdx'dx^ 



(53) 



From Eq. ([24]) and recalling fJ^ ^ i _ 2aip (Eq. p8)) 1. 
the Jordan and Einstein frame metric perturbations are 
related by 



^ ^ ^ + aip , "i/ ^ 'i' ~ aip 

Note how $ ^ vj/ even if l> = ^. 

The Jordan frame metric equation is 



(54) 



[T 



(55) 



where 



(56) 



The scalar field energy-momentum follows from (the min- 
imally coupled portion of) the Jordan frame action Eq. 
([25]) . Note that tp does not have a canonical kinetic term, 
hence the presence of h{(f>) (Eq. (pS)) ). In some theories, 
h might even be zero, such as in f{R). 

The Einstein-Infeld-Hoffmann approach to the prob- 
lem of motion is to compute the gravitational force as a 
surface integral of the momentum flux (Eq. (|3])). The 
momentum flux can be obtained from the spatial com- 
ponents of the pseudo-energy-momentum tensor, which 
is defined by Eq. ([1]): 



SttG 



2 s^^nn^] 



(57) 



Comparing this with Eq. we can see the calcula- 

tion should be similar to the Einstein frame case, except 
that we have an overall conformal factor fl^^ multiplying 
the matter and scalar field energy-momentum, and that 
there are some uniquely Jordan frame contributions to 
tfj^"^ (second line). These differences originate from the 



differences in the metric equation (j55|) from the Einstein 
equation: a non-constant effective 'G", and extra sources 
for Gfi'^ beyond the matter and scalar energy-momentum. 

Phrased in this way, it is not obvious that the integral 
of momentum flux should imply geodesic motion in Jor- 
dan frame. In fact, it does not in general. Relegating the 
details to Appendix [Cl it can be shown that the center 
of mass of an object moves according to: 



MX' = -Af 9i$o + (1 - e)aMdiipo 



(58) 



where $o and (fo are respectively the back- 
ground/environmental (time-time) metric perturbation 
and scalar field, M is the inertial mass and AT* is the 
center of mass coordinate. This is manifestly consistent 
with the Einstein frame equation of motion Eq. ([35]) 
once we recognize that <f>o = 'I'o ^ ct'Po according to Eq. 

([SI. 
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The parameter e is controlled by the relative size of 
the object's gravitational potential and the environmen- 
tal scalar field; see Eqs. (|4T|) & ((42|) . An unscreened 
object, which has e — 1, would move on a geodesic in 
Jordan frame, just as expected for an infinitesimal test 
particle. A screened object, where the chameleon mech- 
anism operates to make e 0, would not move on a 
geodesic even in Jordan frame. In other words, it ap- 
pears as if for screened objects, the gravitational mass 
and the inertial mass are unequal. This can be made 
more explicit if we make a simplifying assumption, that 
both $0 — capo and ipo are sourced mainly by density 
(Eq. dM])), in which case ipo = 2a$o/(l + 2a^), and the 
equation of motion simplifies to 



MX' = M 



1 + 2ea^ 
l + 2a2 



(59) 



which is consistent with its Einstein frame counterpart 
Eq. (|5ip . Thus, in Jordan frame, 



inertial mass = M 
A 

grav. mass 



M- 



(60) 



An unscreened object is subject to a gravitational force 
that is 1 + 20!^ larger than a screened object. For theo- 
ries such as f{R), this is not a small effect: a — 1/V6 
implies 1 20^ = 4/3 (Appendix [Q . This might come 
as a surprise in the Jordan frame, where the metric is 
minimally coupled to matter (Eq. ([25])). Despite con- 
ventional wisdom, specifying the Jordan frame metric 
does not completely specify the motion of an extended 
object. The motion of an object should be ultimately de- 
termined by energy-momentum conservation. Since the 
scalar field contributes energy-momentum, it retains a 
direct infiuence on the motion of the object. This direct 
influence magically cancels only for unscreened objects, 
which thus move on geodesies.® The key physics is the 
(non-relativistic and classical) renormalization of scalar 
charge by the nonlinear chameleon mechanism, as dis- 
cussed in ijlV Al An object's coupling to the background 
scalar field is suppressed if it is screened, and unsup- 
pressed if it is not. ^ 



V. THE PROBLEM OF MOTION IN 
SCALAR-TENSOR THEORIES — SCREENING 
BY STRONG COUPLING A LA VAINSHTEIN 

Another class of scalar-tensor theories where the scalar 
is efficiently screened at short distances is massive grav- 
ity theories. These are not scalar-tensor theories to the 
letter, for at scales of order of the graviton's Compton 
wavelength, and larger, there is no propagating scalar 
degree of freedom. However, at much shorter scales the 
(helicity-zero) longitudinal polarization of the graviton 
behaves like a scalar universally coupled to matter with 
gravitational strength, and the dynamics is accurately 
described by a scalar-tensor theory [l^. What makes 
these theories potentially interesting is that this effec- 
tive scalar comes equipped with important derivative 
self-interactions, which can make it decouple at short 
distances from compact objects, thus recovering agree- 
ment with solar system tests. This is known as the Vain- 
shtein effect, and it shows up explicitly in the DGP model 
— which is a much healthier alternative to the mini- 
mal Fierz-Pauli theory of massive gravitons. However as 
we will see shortly, unlike the chameleon mechanism the 
Vainshtein effect does not lead to 0(1) violations of the 
equivalence principle — in fact it leads to the same vio- 
lations one would have for a purely linear Brans-Dicke 
theory, which are 0(1/0^). 

For simplicity and definiteness, we will analyze the ex- 
ample of DGP. Everything we say can be straightfor- 
wardly extended to the 'galileon', which was argued in 
ref. to be the broadest generalization of (the scalar- 
tensor limit of) DGP that robustly implements the Vain- 
shtein effect without introducing ghost degrees of free- 
dom. The scalar sector of DGP (in Einstein frame) has 
the Lagrangian Q 



-'SAI^iidnf 



M2 

2^(97r)2n7r-h7rT™ ^ 



(61) 



where m is the DGP critical mass scale — the analogue of 
the graviton's mass — and is generally set to the current 
inverse Hubble scale. We are using a dimensionless nor- 
malization for the scalar field tt, so that the Jordan-frame 
metric to which matter fields couple minimally is 



* The cancellation occurs between the effective Newton constant 
Gfl~^ and the uniquely Jordan frame terms Q~^ddQ^ in the 
pseudo-energy-momentum (Eq. i57l ). 

® For a screened object (let's say a galaxy halo) in an unscreened 
environment, it is worth noting that the scalar field ip experiences 
a jump at the thin shell of the object (from a large expectation 
value outside to a small one inside), which means that while 
the Einstein frame metric perturbations are smooth across the 
shell, the Jordan frame metric perturbations jump (Eq. II54I I). 
This resolves what might perhaps be puzzling to some: that the 
screened halo as a whole does not move on the Jordan frame 
geodesic, while infinitesimal particles inside it do. The key ob- 
servation is that they see rather different overall Jordan frame 
metrics (even after accounting for the fact that a particle inside 



hfiL' — hfii^ ~^ 27r rff^i, 



(62) 



That is, in the non-relativistic limit tt is a shift in the 
Newtonian potential: 



$ = $ + TT 



(63) 



the object of course sees additional metric perturbations arising 
from other particles of the object). In Einstein frame, there is no 
mystery: neither the screened halo nor its constituent particles 
experience any scalar force. We thank Wayne Hu for discussions 
on this issue. 
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FIG. 2; Our mathematical trick. 



Comparing Eq. ((6T|) with Eq. ((23)) . we can see that tt 
here is equivalent to ip/^/& there. We adopt the non- 
canonically normalized tt in this section (not to be con- 
fused with the number 3.141...) to ease comparison with 
the DGP literature. In the language of the previous 
section, DGP (and its galileon generalization) has a — 
1 / -v/G, meaning its Jordan frame action has no quadratic 
scalar kinetic term. In this respect, DGP/generalizations 
resemble f{R) or Brans-Dicke theory with ojbd = (Ap- 
pendix [B]), but the resemblance is superficial because 
DGP/gcncralizations have higher derivative interactions 
which allow a completely different screening mechanism. 

A compact source creates a tt profile that scales like 
l/r at large distances (linear regime). However in ap- 
proaching the source the non- hnear termi° of Eq. (pT]) 
starts becoming important, thus changing the behavior 
of TT to ^Jr. This suppresses the force mediated by tt for 
all objects that are inside this non-linear 'halo' surround- 
ing the source. The size of this halo is set by the so-called 
Vainshtein radius, 

where M is the mass of the source, and m is the graviton 
mass scale (about Hubble scale today). 

Despite the qualitative similarities, this kind of non- 
linearity differs significantly from that in the chameleon 
case. The point is that the strong non-linearities tt expe- 
riences close to the source do not renormalize the source's 
total scalar charge. To see this, notice that the source- 



free part of the Lagrangian ((6T|) has a shift-symmetry, 

TT — > TT -|- a . (65) 

As a consequence, tt's equation of motion is the diver- 
gence of the associated Nother's current 

Q^ jf. ^ ^ (66) 

where is a non-linear function of first and second 
derivatives of tt. We therefore have a non-linear Gauss- 
like law for TT, whose source is the trace of the matter 
energy-momentum tensor. For a time-independent solu- 
tion, this equation is readily integrated via Gauss's the- 
orem to yield 

£j ■ da = Q tot (67) 

where S is any surface, and Qtot is the total scalar 
charge enclosed by S, defined as the volume integral 
of (minus) the matter T™" — the total mass for non- 
relativistic sources. So, for instance, if we consider a 
single compact spherical source at the origin and we take 
S" to be a sphere centered at the origin, with radius much 
larger than then on S" tt is in the linear regime, and 
we have 

r ^ 1 GM 

6Mp, <p VTT-dd~M ^ TT ~ r > i?v , 

Js 3 r 

(68) 

regardless of how non-linear tt's dynamics may be in- 
side the Vainshtein region. The mass M is thus also the 



Non-linear in the tt equation of motion, cubic in the action. 



^ Qtot here is proportional to Q in Eq. II43I I. The difference in 
normaUzation is related to the different normalizations for Lp and 

TT. 
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scalar charge, and according to the analysis of mill M 
then is the coupling between the object and very long- 
wavelength external/background tt's. This implies there 
is no chameleon-like non-relativistic renormalization of 
the scalar charge, and no 0(1) violation of the equiv- 
alence principle. The only violations of the equivalence 
principle we may have are associated with how relativistic 
the source's internal structure is, i.e. on how the volume 
integral of T^p** differs from the total ADM mass, which 
is a post-Newtonian 0(l/c^) effect For the benefit of 
the skeptical reader, in Appendix [D] we explicitly check 
this statement along the lines of the previous sections. 

This would be the end of the story, if our universe were 
made up of very isolated compact sources, with typical 
distances between one another much bigger that their 
typical Vainshtein radius. Unfortunately we live precisely 
in the opposite limit. For macroscopic sources the Vain- 
shtein radius is so huge that treating those sources in 
isolation makes no sense. For instance, for the sun we 
have 

i?® - 10^1 cm , (69) 

much bigger than the distance to the closest stars. We 
can repeat the same estimate at all scales, and check 
that nowhere in the universe tt is in the linear regime. 
This makes a comprehensive study of the effects of tt in 
our universe a formidable task, because by definition of a 
non-linear system, we cannot simply superimpose the tt 
fields due to different sources. However here we want to 
focus on possible violations of the equivalence principle, 
and for these we will be able to draw surprisingly general 
conclusions. 

Consider then the more realistic situation of an object 
in the presence of other sources. To be as general as pos- 
sible, let's imagine that the object in question has sizable 
multiple moments as well. The situation is schematically 
depicted in Fig. [2lleft). We would like to compute the 
force acting on the object following the approach of the 
previous sections. In order to approximate the vr field 
created by the other sources (ttq, in the figure) as a lin- 
ear pure-gradient field, we have to draw our surface S 
very close to the object. In fact since tt is deep in the 
non-linear regime, we have two problems 

1. How do we separate ttq from tti — the tt generated 
by the object? The total tt is not simply the sum 
of the two contributions. 

2. What is the effect of the object's multipole mo- 
ments? In a non-linear regime, the multipoles of tt 
are not directly related to the object's multipoles. 
Also, their behavior with r w.r.t. to the monopole's 
is different than in the linear case. 



Notice that since in Eq. H68p we are considering a huge volume 
with sizable non-linearities for tt, we have to make sure that tt's 
contribution to the total ADM mass be negligible. It is straight- 
forward to check that this is indeed the case. 



As to point 1, we are rescued once again by the sym- 
metries of the TT dynamics, which will also make point 2 
uninteresting for our purpose. 

Besides the aforementioned shift-symmetry on vr, the 
source-free part of the Lagrangian (|6ip is invariant under 
a constant shift in the derivative of tt, 

df,TT d^TT + , (70) 

dubbed 'galilean invariance' in ref. [HI, [l^]. In fact 
this is what makes DGP — and suitable generalizations 
thereof — capable of the Vainshtein effect in the presence 
of generic localized sources, and of sustaining interesting 
non-linear cosmological solutions [ll|. Galilean invari- 
ance allows us to add a pure-gradient field to any non- 
linear solution, to get a new non-linear solution. Consider 
therefore the vr field generated by all other sources in the 
absence of our object, and denote this by ttq. Now, if our 
object is much smaller than the typical variation scale of 
ttq, which is set by the distance to the other sources, then 
we can approximate ttq as a constant-gradient field in a 
neighborhood of the object. This means that in such a 
neighborhood, we can freely add ttq to the non-linear so- 
lution vTi we would have for the object in isolation, i.e. in 
the absence of the other sources. The full tt field 

7r = 7ro+7ri (71) 

is a consistent solution to the non-linear problem in the 
neighborhood in question. It is in fact the solution, be- 
cause it has (a) the right source — the object, and (b) 
the right boundary conditions: in moving away from the 
object, the perturbation in S^tt induced by the object 
decays away as l/^/r and we are left with ^^ttq. 

In summary, thanks to Galilean invariance the split- 
ting of Eq. ([7T|l is well defined whenever the object is 
much smaller than its typical distance to nearby sources, 
which is a much milder requirement than linearity. We 
can then draw our sphere S very close to the object, as 
in Fig. mieft), and compute the force acting on it along 
the lines of the previous sections. We do this explicitly 
in Appendix [D] for a spherical object. However for an 
irregular object with sizable multipole moments, we do 
not even know the non-linear solution tti for the object 
in isolation. Nevertheless, here we can use a mathemat- 
ical trick. We have argued that the full solution on S 
is the field the object would have in isolation, tti, plus 
a constant-gradient field, ttq. Nothing prevents us from 
considering a physically different situation, but with the 
same total field on S. This would of course yield the 
same result for the total force acting on the object. A 
very convenient situation to consider is, the object in iso- 
lation with the same background constant-gradient ttq as 
in our case, but extended linearly throughout space, as 
depicted in Fig. [5]Jright). That is, we can forget about 
the external sources, and declare that in fact ttq has a 
constant gradient throughout space. In doing so, we are 
not changing the total tt and its gradient on S — thus we 
are not changing the total force we would compute by 
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integrating over S. We can now change the radius of 
S and bring it outside the Vainshtein region, where it's 
dynamics is hnear, its multipoles decay faster than the 
monopole, and the computation is precisely as for a hnear 
scalar field. By changing the radius of S, we are modify- 
ing the total mass enclosed by S by tt's contribution to it. 
As long as this is negligible, as can be explicitly checked, 
the total force computed this way does not depend on 
the radius of S. 

We therefore conclude that, as far as the coupling of 
an extended object to a background tt field is concerned, 
theories with the Vainshtein effect behave exactly as the 
linear Brans-Dicke theory. The Vainshtein mechanism 
entails equivalence principle violations only of the post- 
Newtonian 0(l/c^) type. 



VI. DISCUSSION — OBSERVATIONAL TESTS 

Let us briefly summarize the main results so far, es- 
pecially those relevant for a discussion of observational 
tests. 

1. In theories where screening of the scalar field occurs 
via strong coupling effects (Vainshtein), such as DGP, 
the coupling of objects to the external scalar fields is un- 
suppressed to good approximation for all objects whose 
mass is not dominated by the gravitational binding en- 
ergy. They therefore fall at the same rate under grav- 
ity. There remains, however, equivalence principle vio- 
lations at the level of <I>, or 0(l/c^) in the language of 
post-Newtonian expansion, just as in Brans-Dicke theory, 
where $ is the gravitational potential depth of the ob- 
ject in question. Order one violation of the equivalence 
principle is visible only when comparing the motion of 
relativistic objects such as black holes (|$| ^ 1) with the 
motion of less compact objects such as stars (j^l ^ 1). 
We will focus on observational tests of the chameleon 
mechanism instead, where the equivalence principle vio- 
lations are more dramatic and readily detectable even in 
non-relativistic / non-compact objects. 

2. The chameleon mechanism implies that the cou- 
pling of screened objects to the external/background 
scalar field is highly suppressed relative to unscreened ob- 
jects. They therefore fall at significantly different rates. 
More precisely, the equation of motion for the center of 
mass of an object is, in Jordan frame (Eq. ([55]) ): 



MX' 



^M^,^o + il-e)aM^^ipo 



(72) 



where X' is the center of mass coordinate, M 
is the inertial mass, $o is the time-time part of 
the external/background metric and ipo is the exter- 
nal/background scalar field (through which the object 
moves). Here we have adopted a convenient parame- 
terization of the relevant conformal factor ^^^(9?) in the 
scalar-tensor action (Eqs. & 



(73) 



where a is a constant. Eq. (|73p can be thought of as 
a small field expansion of e.g. exp[— 2a(^]. The quantity 
e represents the degree of screening. Unscreened objects 
have e = 1, while screened objects have e < 1. The 
screening parameter e is controlled by (Eqs. ((4T|) & ((42|) ): 



2a{GM/rc) 
e~ 1 ff 



if 



2a{GM/r,] 



< 1 (screened) 



V5* 



2a{GM/ra) 



>1 (unscreened) (74) 



where Tj, is the size of the object in question, and ip^ is the 
value of the external scalar field. Screening is therefore 
controlled by the ratio of the exterior scalar field and the 
gravitational potential of the object. Unless the values of 
the two are finely tuned, one generally expects e ^ 1 if 
screening occurs at all. Eq. ([7^ tells us only unscreened 
objects, those with a sufficiently shallow gravitational po- 
tential, have equal inertial mass and gravitational mass 
in Jordan frame. 

3. The precise ratio of inertial mass to gravitational 
mass predicted by Eq. ([7^ depends on how the external 
scalar field (po is sourced. There are two useful limits to 
consider. One is when the external scalar field is sourced 
primarily by density (Eq. (|C17[) ) in which case the equa- 
tion of motion is well-approximated by 



MX' = -M 



1 + 2ea^ 
l + 2a^ 



(75) 



n^{(p) = l-2ap , I 1 



which implies the unscreened and screened gravitational 
mass are related by 

unscreened g. mass = (1 + 2q!^) x screened g. mass(76) 

The other limiting case is when the external scalar field 
profile is determined primarily by the potential, in which 
case ipo is Yukawa suppressed on sufficiently large scales 
(scales larger than the Compton wavelength of the scalar 
field ~ 1/m^). If ipo is Yukawa suppressed, the sec- 
ond term on the right of Eq. ([7^ is small regardless 
of whether the object is screened or not. In this regime, 
one cannot observe equivalence principle violations by 
comparing motions of different objects. As we will dis- 
cuss below, it is therefore important to avoid Yukawa 
suppression when performing observational tests. 

4. Observationally, the single most important parame- 
ter besides in a chameleon theory is the quantity 
(see Eq. ([731 ). For a given object, the relevant (p* is the 
scalar field value in its immediate surrounding. If this 
is sufficiently small, the object will be screened accord- 
ing to Eq. ([7^ . In other words, if the object is sitting 
in some high density environment which is itself screened 
(i.e. small scalar field value), it is likely that the object is 
also screened. We will refer to this as blanket screening. 
In an environment where blanket screening of all objects 
applies, violations of equivalence principle will be hard 
to detect. 

Conversely, if the object is sitting in an environment 
where the density is at the cosmic mean or even lower 
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(voids), the chances are better that the object is un- 
screened because the relevant is larger. In other 
words, to find unscreened objects, it is useful to avoid 
overdense regions and to focus on objects with a shallow 
gravitational potential such that (Eq. (|74p ): 



GM ^ ip^ 



(77) 



where is taken to be the cosmic mean value. We 
already know ip:t/{2a) < 10^® from requiring that the 
Milky Way be screened to avoid gross deviations from 
GR in solar system experiments [22| (see also [i^ for 
a recent revision of the Milky Way's mass). This 
leaves us an interesting window for testing theories like 
f{R)- The smallest galaxies of interest have GM/tc ^ 
10~^(w/30 km/s)^. Therefore, if for our modified 

gravity theory satisfies 
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-8 < ^ < IQ- 
^ 2a ^ 



(observationally interesting) (78) 



massive galaxies at the high end would be screened while 
less massive galaxies at the low end would not, offer- 
ing many opportunities for observing equivalence princi- 
ple violations, as we will see. The parameter space for 
chameleon mechanism is illustrated in Fig. [31 
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FIG. 3: Schematic experimental bounds in the chameleon pa- 
rameter space. The shaded region is already excluded, by 
demanding that the Milky way be screened. The dashed line 
is the improvement we propose. Further improvement is pos- 
sible using Lyman-alpha clouds. For a <C 1 the chameleon 
is very weakly coupled to matter to begin with, and clearly 
this relaxes the bound on the other parameter. Curvature 
invariant theories such as f{R) have a = l/\/6, though this 
value is not protected by symmetry. The generic expectation 
is Q ~ Oil). 

5. To make our discussion more concrete, it is useful 
to keep in mind a particular class of examples, the so 
called f{R) theories. As shown in Appendix |b1 f{R) 
corresponds to a specific value for a: 



a = l/Ve for f{R) 



and so the unscreened gravitational force is 4/3 higher 
than the screened one (Eq. ([75]) ). This value for a arises 
simply from demanding that ip be an auxiliary field in 
Jordan frame (i.e. h{if) = in Eq. (HI])). It is worth 
stressing however that from a theoretical viewpoint f{R) 
theories are not favored over more generic scalar-theories. 
For instance the property that Lp has no kinetic energy 
in Jordan frame is not protected by any symmetry, thus 
one generically expects that upon quantum corrections 
an /(i?) theory will become a scalar-tensor theory not be- 
longing to the f{R) class. In fact, the resulting a could be 
even larger than the 1 / On the other hand, the struc- 
ture of the TT Lagrangian in DGP and in galileon theories 
is protected by a symmetry — Galilean symmetry — which 
forbids the quantum generation of Lagrangian terms with 
fewer than two derivatives per field |9]. 

The most important parameter in an f{R) theory is 
fii = df/ dR evaluated at the cosmic mean, and it can be 
shown 

Jr = ~2ap, . (80) 

Appendix |B] also shows that common parameterizations 
of f{R) theories connect (p^, and the Compton wavelength 
l/m^p of the scalar field: 



-H-^ (81) 



where H^^ is the Hubble scale. The fifth force is Yukawa 
suppressed between a pair of objects separated by scales 
larger than This places an important restriction 

on the range of observational tests available for these the- 
ories. One should keep in mind though that Eq. ((8T|) can 
be evaded by considering more general potentials than is 
commonly discussed (see Appendix [Bjl . 

Observational tests. Below we list a series of possi- 
ble observational tests i5Q|]- They can be viewed in one 
of two ways: if they turn up positive, we will have a mea- 
surement of ip^:/a and a within the chameleon framework 
and a refutation of screening by strong coupling; if they 
turn up negative, we can interpret the results either as an 
upper limit on (p^/a and a for the chameleon mechanism 
or as indirect support for the strong coupling picture (or 
as yet another confirmation of GR of course) . Tests 1 and 
2 focus on the bulk motion of individual galaxies while 
tests 3, 4 and 5 focus on the internal motion of their con- 
stituents. Fig. m illustrates these various possibilities. 
Further tests are discussed at the end. 

Test 1. Eq. (|75p tells us that unscreened galax- 
ies (whose gravitational potential is sufhciently shallow 
that Eq. ([77)1 is satisfied) accelerate faster compared to 
screened galaxies by a factor of 1 -|- 2a^, which is 4/3 
for f{R) theories. We therefore expect small galaxies to 
move faster than large ones in general. One must be 
careful not to confuse this with dynamical friction which 
could cause a similar effect (sometimes referred to as ve- 
locity bias; see [5ll.[5^). 

A common way to probe the peculiar motion of galax- 
(79) ies is via redshift space distortions. On large (linear) 
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scales, where dynamical friction is unimportant, one 
could use the anisotropy of the redshift space correla- 
tion function to measure /3 — b^^dlnD/dlna where b 
is the linear galaxy bias, D is the linear growth factor 
and a is the scale factor [s^ [s^]- The new twist intro- 
duced by the chameleon mechanism is that effectively 
d\n D/dlna, which arises from the ratio between pecu- 
liar motion and mass overdensity, could be different for 
different size galaxies. This means if one takes the ratio 
/3i//32 between small and large galaxies, and multiplies it 
by the appropriate ratio of the real space power spectrum 
P to remove the dependence on galaxy bias, GR would 
return unity while the chameleon mechanism would not 
i.e. 



7^ 1 for chameleon 



(82) 



where 1 denotes the unscreened (small) galaxies and 2 
denotes the screened (large) galaxies. The precise pre- 
diction for chameleon theory is a bit complex and will be 
presented elsewhere. This is a worthwhile test to carry 
out — if anything, it serves as a consistency test for GR — 
but as a test of f{R) theories, one must take care to avoid 
Yukawa suppression. For those classes of f{R) theories 
that obey Eq. ([5T|) . the Compton wavelength is fairly 
short given current upper limits on ipi,/a (Eq. (j78[) ). In 
that case, Yukawa suppression implies the relative mo- 
tion between pairs of galaxies separated on linear scales 
is likely insensitive to the scalar force. We emphasize, 
however, that the relation between Compton wavelength 
and (f^/a in Eq. ([8T|) is not general and can be relaxed. 

Alternatively, one could focus on the relative pecu- 
liar motion between galaxies on small scales, but then 
one needs to contend with the possibility of both blan- 
ket screening and dynamical friction. One way around 
these problems is to adopt the method of ^bE] which al- 
lows one to study the pairwise velocity dispersion as a 
function of environment. Voids would be the ideal envi- 
ronment which bypasses both problems, and allows one 
in principle to observe the systematically different pecu- 
liar velocities between screened and unscreened galaxies. 
However, because the observed peculiar motion is always 
galaxy weighted, one must devise a way to disentangle 
apparent velocity difference due to different (nonlinear) 
galaxy biases from genuine velocity difference due to the 
chameleon mechanism. 

In summary, the cleanest galaxy bulk motion test is 
probably to examine redshift distortions on large (linear) 
scales where there are no complications from dynamical 
friction or nonlinear galaxy bias. In this case, a neces- 
sary requirement for order unity observable effects is that 



the scalar Compton wavelength (at cosmic mean density) 
cannot be too small, i.e. it must extend to linear scales. 

Test 2. In environments such as voids, where blanket 
screening of all galaxies docs not occur, the systematic 
difference in the motion of small versus large galaxies im- 
plies their density distributions would differ from stan- 
dard expectations. For instance, galaxies are expected 
to stream away from voids. The fact that small galax- 
ies stream at a faster velocity means voids defined by 
them would appear larger than expectations based on the 
large galaxies. This prediction of the chameleon mecha- 
nism is intriguing because as emphasized by Peebles (sGj 
the observed voids defined by the small galaxies appear 
too large compared to predictions from ACDM model, 
whereas the observed void sizes defined by larger galax- 
ies appear consistent with the standard model [57 . 58| . 
For the chameleon mechanism to provide a viable solu- 
tion to the void problem, one must consider models that 
evade the Compton wavelength condition of Eq. (|5T|) . 

More generally, because screened and unscreened 
galaxies move differently, their clustering bias will evolve 
differently, giving potentially observable effects which 
will be explored elsewhere (see e.g.jsoj). 

Test 3. Perhaps the cleanest way to test the 
chameleon mechanism is to examine the internal motions 
of objects inside an unscreened galaxy i.e. whose poten- 
tial is sufhciently shallow that GM/rc < ip*/{2a) (Eq. 
((77)) ). The sun has 



GM 



2 X 10"^ 



(83) 



and main sequence stars are expected to have a sim- 
ilar gravitational potential [60j . Current constraints 
{^*/{2a) < 10^6; Eq. jTH])) therefore likely imply the 
majority of stars inside an unscreened galaxy are them- 
selves screened^ because the stars have a fairly deep po- 
tential. On the other hand, diffuse gas inside the same 
galaxy likely remains unscreened because the gravita- 
tional potential at a gas cloud should be dominated by 
the potential of the galaxy rather than that of the gas 
cloud itself. For instance HI gas clouds typically have an 
atomic number density of ~ 10 cm^"^ and a size of ~ 10 pc 
[6l| giving a tiny self-potential of GM/rc ^ 10^^^. 
Therefore, at the same radius from galactic center, the 
HI gas clouds should move with an acceleration that is 
larger than that of the stars by a factor of 1 -|- 2a^, which 
is 4/3 for f{R). Mass estimates constructed from these 
two tracers will differ by the same factor By focusing 
on internal motions, this test avoids issues to do with 



This can be obtained from looking at Fourier modes that point 
transverse to the hne of sight. 



One might also worry that a galactic disk could have a deeper 
gravitational potential than the halo and potentially blanket 
screen all objects in it. In practice, the gravitational potential of 
a disk is comparable to that of the halo [63 |. 
When comparing HI rotational velocity with stellar rotational 
velocity from a disk galaxy, care should be taken to account for 
a possible asymmetric drift 63j. A proper dynamical modeling 
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FIG. 4: A schematic illustration of observational tests. $acif 
and 'I'cnv represent the gravitational potentials of an object 
and its environment. They can be thought of as ~ {v/cf', 
where v is the internal velocity. The value (^t/(2a) delineates 
screening or lack thereof — objects/environments with a po- 
tential deeper than <p*/(2a) are screened (shaded), and those 
with a shallower potential are unscreened (unshaded). Cur- 
rent constraints tell us (pt/(2a) has to be less than ~ 10~®. 
There are many comparison tests one can make. For instance, 
an unscreened diffuse gas cloud residing in a dwarf galaxy 
(A), versus a screened star residing in the same galaxy (B) — 
A falls faster than B. This situation can be replicated on a 
larger scale e.g. a dwarf galaxy in the fields/voids (A) ver- 
sus a massive galaxy in the same fields/voids (B). Another 
example: a dwarf galaxy out in the fields/voids (A), versus a 
dwarf galaxy residing in a group or cluster (D) — D is blanket 
screened by its environment and would exhibit no equivalence 
principle violations in its internal motions of gas clouds and 
stars, while A would have observable violations. On the other 
hand, a massive galaxy would have no such (internal) equiva- 
lence principle violations whether it be in the fields/voids (B) 
or in a group/cluster (C). 



Yukawa suppression that plagues the bulk motion tests 
of some /(i?) models i.e. the scalar Compton wavelength 
probably exceeds the size of a galaxy (Eq. (|8T|) ) . To max- 
imize the chance that the galaxy studied is unscreened, 
one should look for the smallest galaxies, preferably in 
voids or at least in the fields A small galaxy with 



should unveil systematic difference in the mass estimated from 
HI and that from stars, should the chameleon mechanism be at 
work. 

Screened galaxies are not interesting for this test because both 
the stars and the HI gas would be blanket screened. 



internal v ~ 30km/s would provide us either a positive 
detection of the chameleon mechanism, or an upper limit 
on (p^/{2a) of about 10~®, for an 0(1) a like in f{R). 

Test 4. In addition to estimating masses of galax- 
ies from their internal dynamics, one could also esti- 
mate their masses using gravitational lensing. Photons 
should behave as unscreened particles and move on null 
geodesies in the Jordan frame. They therefore see the 
($ + ^')/2 potential. Unscreened non-relativistic objects 
(such as HI gas clouds) on the other hand should move ac- 
cording to the $ potential, while screened non-relativistic 
objects (such as stars) should move according to Eq. ([7^ 
with e '--^ 0. Eq. (|C6p together with the assumption of 
negligible Yukawa suppression then tells us that the mass 
estimate from photons should equal the mass estimate 
from stars, and both should be smaller than the mass 
estimate from HI by a factor of 1 + 2a^ . Carrying out 
this test might be a bit challenging as one would ideally 
like to measure the lensing mass for the smallest galaxies 
possible to avoid screening of the whole galaxy. Stack- 
ing small galaxies and performing a weak lensing shear 
measurement is likely the way to go. 

Incidentally, returning briefly to the subject of bulk 
motion: for the same reasons articulated above, weak 
lensing measurements on large scales should yield con- 
sistent results when compared against redshift distortion 
measurements of screened (large) galaxies, but inconsis- 
tent results when compared against redshift distortions 
of unscreened (small) galaxies This adds a possible 
new twist to existing tests of this sort proposed in the 
literature [63 |. 

Test 5. Kesden and Kamionkowski [6^ recently pro- 
posed an interesting test of the equivalence principle: if 
dark matter particles and stars fall at different rates un- 
der gravity, the tidal streams from infalling galactic satel- 
lites would be asymmetric. To this the chameleon mech- 
anism adds a new twist: only unscreened (small) galax- 
ies will exhibit this asymmetry; screened (large) galaxies 
will not, because both dark matter and stars are blanket 
screened. The current limit of (y9,/(2a) < 10^^ guaran- 
tees that the Milky Way does not exhibit the Kesden- 
Kamionkowski effect. To see it, one would have to ex- 
amine smaller galaxies, preferably out in the fields or 
voids. It is also important to emphasize the difference in 
the origin of equivalence principle violation in this paper 
from that in most of the cosmology literature, including 
[65| . In most of the literature, equivalence principle vi- 
olation is there from the beginning, in the sense that it 
exists at the level of the microscopic action i.e. elemen- 
tary baryons and dark matter particles are coupled to 
the scalar differently (see discussion in ^|l|. In this pa- 
per, there is no explicit equivalence principle violation at 
the level of the microscopic action — it gets turned on by 
the chameleon mechanism and differentiates not so much 



Modulo Yukawa suppression, as usual. 
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between baryons and dark matter, but between objects 
with a deep gravitational potential and those with a shal- 
low potential. 

The tidal stream test brings to mind an interesting 
question: if stars and dark matter inside an unscreened 
galaxy fall at different rates under some external gravi- 
tational field, shouldn't the galaxy fall apart in the sense 
that the stars and dark matter separate (even in the ab- 
sence of significant tidal forces)? This should happen 
at some level, but should not be a large effect as long as 
internal gravitational forces dominate over external ones. 

Perhaps the ultimate test objects are low column den- 
sity Lyman-alpha clouds in the intergalactic medium. 
They have very shallow gravitational potentials and can 
therefore be used to probe a much lower Lp^/{2a) than 
are afforded by dwarf galaxies. Comparing the motions of 
clouds [6^ . [67j and galaxies in similar environments could 
unveil the chameleon mechanism at work or put severe 
limits beyond those depicted in Fig. [S] This is probably 
a good place to end this discussion. Many of the ideas 
presented above need to be worked out in detail, and 
undoubtedly there are many other possible observational 
tests. We hope to pursue these in the future. 
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APPENDIX A: THE EFFECT OF HIGHER 
MULTIPOLES 

We want to show that as long as the external fields 
$Oj^o are approximately linear on the spherical sur- 
face S*, higher multipole moments of the object have no 
effect — as expected, tidal effects are there only if the ex- 
ternal field gradients vary appreciably over the object's 
volume. We will do this explicitly for GR. The gener- 
alization to the scalar-tensor case is straightforward in 
Einstein frame. Our discussion here thus makes clear the 
surface S can in fact be chosen quite close to the object; 
it is not necessary to choose a large radius r so that only 
the monopole survives. As we will see, even the spherical 
shape of the surface is unnecessary. So, S can just be the 
object's surface. 

Let us not assume that the object's multipole moments 
are negligible. We have to dot Eq. ^ with dSj, and in- 
tegrate over a surface S surrounding the object. We need 
not assume this surface is spherical. The first simplifica- 



tion is that (Eq. ([371) 1 

$ = * , (Al) 

which means that the term proportional to G^^\^ van- 
ishes, and that many terms are equal. Also 

V^* = , (A2) 

as long as the environmental p right at S is negligible 
(which it must be if we are approximating $0 = 5*0 as 
linear). We are left with 

G(2)iJ = 4$ d,dj<^ + 29,$ dj<^ - iS^jidk'^f . (A3) 

The second simplification comes from the fact that we 
have to integrate this in dSj, over a closed surface. This 
means that we can add a curl for free. We use this free- 
dom to eliminate the second-derivative term, thus getting 

G(2),^j ^ _2a,$ dj^ + % {dk^f . (A4) 

The total force thus is 

P = - i [2E{E ■ fi) - fiE'^] dS , ( A5) 
Js 

where we defined the 'electric field' 

£; = V$ , (A6) 

and h is the outgoing normal. We now decompose E: 

E = Eo + Ei, (A7) 

where the first term is the environmental background, 

Eq ~ const , (A8) 

and the second is due to the object itself. Expanding 
Eq. (|X5|) we get 

- (SttG) P ^ j> [^EoiEo ■ n) - nE^]dS 

2 ^ [Eo{Ei ■ h) + Ei{Eq ■ h) - uEq ■ Ei]dS 
+ j [2Ei{Ei ■ h) - hEl]dS . (A9) 

The first line vanishes for constant Eq because 

j> h^dS = , (AlO) 

for any closed surface. To show that the last line also 
vanishes in general, we make use of a vector form of the 
Kirchoff integral (see e.g. ref. [S^): if E is an electrostatic 
field with no source on S, and G{x — x') is a Green's 
function, 

V^G{x-x') = S^{x-x') , (All) 
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then 



where 



2E{VG-n)-hE-VG+hx(VGxE) dS = , (A12) 



whenever the point x' is inside the volume bounded by 
S. To be totally clear, _E is a function of x, over which 
we are integrating. G instead is a function both of x and 
of x'. The gradients are taken w.r.t. x. If we replace 
E El, we multiply by the object's density pi{x'), and 
integrate in x' over the object's volume, the l.h.s. reduces 
precisely to the last line of eq. (jA9p . Physically, its van- 
ishing means that the object does not exert any net force 
on itself — as expected. 

We are thus left with the second line of Eq. (|A9[) — the 
mixed terms. Factoring Eq out of the integral, we have 
to compute the tensor integral 



i E{ F dS= <L dS , (A13) 

Js Js 

and then take suitable contractions: 
47rGr 

However, although not manifestly, is a symmetric ten- 
sor: 



d,^if^dS= / d,d,<^id^x 



(A15) 



This means that only the first term in brackets in 
Eq. (|A14|) survives, and we have 



(A16) 



The flux of V(f>i through S defines the total mass en- 
closed by S, and we finally get 



P = -MV4>o 



(A17) 



independently of higher multipole moments of the ob- 
ject's mass distribution. 



4a2 



(B3) 



Phenomenologically allowed Brans-Dicke theory has a 
large wbd which corresponds to a small a. With the ad- 
dition of a suitable potential, the chameleon mechanism 
makes even small wbd theories observationally viable. In- 
deed, f{R) theories correspond to wbd = 0, or a = l/\/6. 
More concretely, consider a scalar-tensor theory with the 
following action parameters (Eq. (^5)1 ): 



1 



df 



0C 



df 



f 



(B4) 



where / is some function of a field cj)c, which is related 
to our familiar (p through 



n'^ = exp [-2aLp] = 1 + 



_df_ 

d(j>c 



(B5) 



with a set to 1/V6. Chiba showed that this is exactly 
equivalent to an f{R) theory with 

s^jd^xv-9^[R+fm 

) • (B6) 



This can be easily seen by integrating out (/)c, or equiva- 
lently (p, which is a pure auxiliary field. Eq. (jBSp makes 
clear why in f{R) theory, df /dR, which is simply df /d(j)c 
for ipc — R, plays the role of a scalar field — for practical 
purposes, fi^ ~ 1 — 2aip, and therefore 



df_ 
dR 



—2a If 



(B7) 



which we also refer to as —2q in Appendix [Cl (Eq . (|C1|) '). 
The Compton wavelength can be derived from 



APPENDIX B: MAPPING TO BRANS-DICKE 

AND f{R) 

It is useful to relate the Jordan frame scalar-tensor ac- 
tion of Eq. (pS)) to certain classes of scalar-tensor theories 
in the literature. 

Brans-Dicke theory, as it is usually written, corre- 
sponds to choosing the potential = 0, and defining 



i>B-D = —^'^{(p) = — exp [-2aip] 



(Bl) 



Eq. (pS)) then becomes 
1 



S = / d-'xy/^- 



167r 

+ d X Cm{lp 



BD 



R 



t^BD 



tfBD 



V^^bdV^^bd 



(B2) 



(B8) 



Note our somewhat unusual normalizations for V and 
(f in Eq. (|^5|) : the actual potential is Mp^V and the 
canonically normalized scalar field is Mpiip. These extra 
factors of Mpi nicely cancel in the definition for the scalar 
mass m^p. For certain potentials, such as a power law, one 
can approximate 



IdV 

p dip 



(B9) 



Assuming ip is sitting in equilibrium, i.e. the right hand 
side of Eq. p9|) vanishes, we obtain 



TO,, 



^p* 



(BIO) 
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where we have assumed the local density p is at the cos- 
mic mean, and is the corresponding scalar field value. 
The Compton wavelength in other environments can be 
computed in an analogous fashion. 

We should emphasize that Eqs. (jB9p and thus (|B10[) 
are by no means forced upon us. A common parameteri- 
zation of the /(i?) potential does have this form f23|, but 
it is not necessary. For instance, an exponential potential 
of the form V oc expP/(p can evade these restrictions, if 
a suitable choice of the coefficient /? is made. 



APPENDIX C: JORDAN FRAME DERIVATION 
— SCREENING BY CHAMELEON 



Here, we supply details that are skipped in i jlVBI for 
our Jordan frame derivation of the equation of motion. 
It will prove convenient to define a parameter q: 



D,'^{if) — l — 2q{(p) i.e. q — aip 



(CI) 



where ^"^{(p) is the conformal transformation between 
Einstein and Jordan frames (see Eqs. p4)) & ([281)). We 
assume |g| <C 1. 

In addition to the Jordan frame metric equation ()55|1 . 
we need the scalar field equation: 

„^+lfiv<.^V>-^ + if!«^0, (C2) 

2o^p dtp 2 dip 

Both the metric equation and the scalar field equation 
look considerably more complicated than their Einstein 
frame counterparts (Eqs. ((30|) & ((32l) ). The differences 
will be crucial in reconciling the Einstein frame and the 
Jordan frame viewpoints on the problem of motion. 

Linearizing the metric perturbations, but not the 
scalar field perturbation, we obtain the following metric 
equations: 

V^* = n-"^ (^^TiGp - i/iVo(^VV + \h\I,p\I'p 

dod^^ = -^n-^hv°pVip>- 
(a,a,-i%v2)(vi'-$) = 

{hv.v^p - ^6i'hVky^v''ip) 

{8nGp + hV^ipV^ip + in^v) 

+3n-^Dn^ (C3) 

where the last equation comes from the space-time trace. 
Combining this trace equation with the scalar field equa- 
tion (Eq. (jC2[) ). we obtain: 

dip dip dip 



This is nothing other than the conformally transformed 
version of the scalar field equation in Einstein frame ( Eq. 
([32]) ). Ignoring time derivatives, metric perturbations 
and using Eqs. ^Si) & (ESI), we obtain 



(C5) 



V^cp = — +a8TrGp . 
dp 



This scalar field equation has exactly the same form as 
its Einstein frame counterpart Eq. (|39|) . and therefore 
the chameleon mechanism works in the same way. Let 
us stress that the important chameleon nonlinearity is 
retained in Eq. (|C5[) . namely that the potential V is not 
linearized. For instance, it could have the form of l/p" 
where n is positive (see Fig. [Ij. 

The metric equations Eq. (jC3|) can be simplified by 
invoking the same set of assumptions we used in the Ein- 
stein frame i.e. time derivatives are small, </? <C 1, and V 
can be ignored both inside the object of interest (because 
Gp ^ V) and outside at the radius of interest r: 

+ $)/2 = inGn'^p ~ 47rGp 
V'^i^ - $) = n-^V^n'^ ~ -2V2g . 

This is the analog of Eq. ((37| in Einstein frame. 

Employing the same decomposition of fields into con- 
tributions from environment (subscript 0) and object 
(subscript 1) as before, 

$ = $0 + *i(r) , * = I'd + *i(r) , (/? = </Jo + <^i(7-)(C7) 
we obtain the following solutions to Eqs. (jCSP & (|C6|) : 



«'o-«'o(0) + a,*o(OK 



(C6) 



*1 = 


GM 


+ 91 




r 


*1 = 


GM 






r 




2GM 


Pi = 


~e a 



(C8) 



where qi is defined below (Eq. (jClip ). These solutions 
are obtained under the same set of assumptions enumer- 
ated at the end of i )IV Al namely we can find a radius r 
outside the object which is: 1. small enough such that 
second derivatives of the environmental fields can be ig- 
nored; 2. large enough such that the monopole contribu- 
tion from the object itself dominates. The solution for ip 
is non-trivial: see discussion in i jIV Al on the chameleon 
mechanism. The quantity e signifies whether chameleon 
screening occurs: it equals unity if there is no screening, 
and becomes small if there is screening (see Eqs. (pi & 



Note that we could have easily obtained these expres- 
sions by applying the small conformal transformation Eq. 
([51)1 to the corresponding Einstein frame results. 

The next step is to use these solutions to compute the 
gravitational force by way of a surface integral at radius 
r (Eq. (HI)). Here, we need to be careful in defining the 
relevant pseudo energy-momentum tensor. The Jordan 
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frame metric equation (Eq. ([551)) can be rewritten in the 
form of Eq. ^ with 



m V 



SttG ^ 



8ttG 



(C9) 



As usual, we ignore the exterior matter energy- 
momentum at radius r. We have aheady computed the 
relevant surface integral involving G*^^-' (Eq. (fTT|) ). We 
need to redo the surface integral over the scalar field 
energy-momentum because unlike in Einstein frame, our 
scalar field here is non-canonical. We also need to con- 
sider a scalar field contribution that is unique to Jordan 
frame (second line of Eq. (IC9P ). 

Eqs. dSni) together with ^ tell us that 



-6[d,qd^q~S,^dkqd''q] } 



(CIO) 



where we have ignored time derivatives, and used fl^ = 
l — 2q with |<7| <^ 1. As expressed in Eq. (^5]) . we can ap- 
proximate dq/ dip as constant, which means we can Taylor 
expand q{'~p) to first order: 



(7 = (7o + 91 where 



qi=aipi (Cll) 



where is some typical scalar field value exterior to the 
object, and (po and ipi are as defined in Eq. (jC8|) . 

Applying the same arguments that go into obtaining 
the Einstein frame result Eq. we find 



OiipQ— \-QOiqa — 

or or 



(C12) 

Lastly, the contribution from the uniquely Jordan 
frame terms (last line of Eq. (jC9ll ) works out to be 



Q.-^ r- 

dSi--^ [ViV^f^^ - 5^^un'^] = — X 



[1 

3G' 



-4.diqo— + Wiqo— o^qo-^ o^<^a — 

or or or or 



(.C13) 



Combining together all these contributions from Eqs. 
dm), (ICT2l) & (IUl3l) . we find 



where for the last equality we have used the fact that 
$1 - ^-1 = 2qi from Eq. (jUS)) . 

Substituting the solutions for $1, and (pi from Eq. 
(jCSp . and equating Pi with MX* where X' is the center 
of mass coordinate (Eq. (fT5|) ). we obtain 

MX' ^ -Md.i^a - go) - eaMdiipo . (C15) 

This result is consistent with the Einstein frame equation 
of motion Eq. (j48|) once we recognize that <i>o = $0 — Qo 
(see Eq. (f54| ). In fact, we can simplify this further by 
using Eq. (jClip : 

MX* = -Mdi<Po + (1 - e)aMdip>o (C16) 



which reproduces Eq. ((58)) . 

We can simplify this further if we make an additional 
assumption about how the environmental fields $0 and 
Po are sourced, which does not necessarily hold in gen- 
eral. That both $0 — qo and po are sourced primarily by 
the density: 

V2($o - go) = 47rGp , VVo = a87rGp (CI7) 

where p is the environmental density. These are the 
analogs of Eq. (|49p in the Einstein frame. With this 
assumption, we can see that 



which then implies Eqs. (|59p and ([60|l . 



(C18) 



APPENDIX D: TOTAL FORCE IN THE 
VAINSHTEIN CASE 

Consider the effective 4D description of DGP. It is a 
scalar-tensor theory, where the scalar Lagrangian in Ein- 
stein frame is given by Eq. (|6ip . Assuming that we have 
a spherical object, we want to compute the total force 
acting on it as 



(Dl) 



where is a sphere surrounding the object, possibly deep 
in the non-linear regime for tt. The effective stress-tensor 
ti^ receives three kinds of contributions: 

1. A gravitational one, bilinear in $ and 5*, propor- 
tional to G(2),J; 



" 2G 



d 



$1 



\<ll) - ^^P:>-Q^ 



r 

2G 



-9.90^0 + 1*1 -f$i) 

' ^ o , 9($i + *i) „ dpi 
dti^o - qo) K d,po 



dr 



dr 



(CM) 



2. A free scalar-like one, schematically of the form 
81:87: ^ coming from the quadratic part of eq. (|6ip : 

3. A cubic one, of the form 8it 8t: 8^1:, coming from 
the cubic part of eq. (|6T|l. 



The first two contributions are already dealt with in i jlll 
^IV Al and we need not compute them again here. We 
will focus instead on the third contribution, which is the 
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intrinsically new feature introduced by the Vainshtein ef- 
fect. 

The cubic part of the stress-energy tensor is 



2M 



PI 



where R is the radius of S. To this, we have to add 
the GR contribution \12\ . as well as the scalar contribu- 
tion coming from the quadratic part of the Lagrangian — 
eq. (Uni), corrected by a factor 6 to make up for the non- 
canonical normalization of tt: 



■r]^,^da7Td°'{dTry 



(D2) 



P = -MV$o - SttMIi i?2 



/2 



37ri 



4 Tr[ 
m? R 



Vtto . (Dll) 



This comes from varying the cubic part of the tt action 
([6T|) w.r.t. to the metric. (The variation of Dtt is straight- 



forward upon rewriting Dtt = --^=df_i{^/—gg^'^ d^ir).) We 
want to compute 



We now have to relate tt'^ to the mass of the source 
and to R. As discussed in ^Vl the field equation deriving 
from (1611) can be rewritten as a Gauss law, 



J s 



(D3) 



We split the total tt field as the sum of an environmental 
one and of that created by the object 



J-dS , 



(D12) 



(D4) 



where for a spherically symmetric 7ri(r) the current J is 



As discussed in ^ this splitting is well defined even in 
the non-linear regime, as long as the external field ttq can 
be approximated as a constant gradient field: 



Mil 



6tt[ 



(D13) 



Vttq ~ const. 



(D5) 



Also, assuming spherical symmetry of the source, we re- 
strict to a radial tti : 



This is precisely the combination appearing in eq. (jPlip . 
We thus see that the total force acting on the object is 
simply 



TTi = 7ri(r) 



(D6) 



P = -AfV$o - MVtto 



(D14) 



Given ()D5[) . second derivatives in (|D2p have to act on tti. 
Also, upon integrating over the sphere the only surviving 
terms in (|D2p are those involving one ttq field and two 
TTi 's — all the others integrate to zero because of spherical 
symmetry. We thus have 



independently of R, that is, independently of whether 
the sphere S is in the linear or non-linear region. This 
equation is in Einstein frame. Transforming to Jordan 
frame gives simply: 



21/2; 



P = -MV$i 







(D15) 



2 afc7ro(aj7ri djdkTTi - Si-j 9/7ri dkdini) 



+ {i^ j) 



(D7) 



Using 



tt' 

didjTTi = {6ij - fifj) -i -I- f,fj tt" 

(primes denote derivatives w.r.t. r) we finally get 

1 r. 



(D8) 
(D9) 



dS- 



-|-9j7ro67ri fifj 

5— -iUTT-^ djTTo , 



(DIO) 



APPENDIX E: RESULTS IN AN FRW UNIVERSE 

Here, we wish to provide extensions of our results to an 
expanding universe. We will do this in Einstein frame. 
Jordan frame results can be obtained by making a con- 
formal transformation (Eqs. ([24]), (IM]) ). 

The metric is 



ds^ = a{rif -{1 + 2^)dr]^ + (1 - 2^)Sijdx'dx^ 

'(El) 

We will use ' for derivative with respect to conformal 
time 77. 

The components of the Einstein tensor up to first order 
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in metric perturbations are 



Go 



Go — 



G/ 
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\ a J 
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a 



a 




(E2) 



The affine connection components are, up to first or- 
der: 



^ 00 


a' 
a 




pO _ 
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- 2* - 2$) - vj/' 






r'oo 




r'oj - % 
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Sjkdi-^ ■ 


(E3) 



Using a ~ to denote the spatially homogeneous compo- 
nents, we have the following zeroth-order Einstein equa- 
tions 



3— r = SvrG 
a* 



^ - -— r = 47rG 
a3 2 a-i 



,(E4) 



The somewhat unusual appearance of factors of M|[ 



arises from our choice of action (Eq. which makes (p 

dimensionless. Working within the same approximation 
as before, we can ignore the kinetic term for ip. The po- 
tential V consists of a constant piece (cosmological con- 
stant) plus a ip dependent piece. The latter is constrained 
by the last expression of Eq. , which tells us if (p is 
sitting at equilibrium, then 5V ^ ap(p/Mp^ ^ p/Mp^. 
This means the non-constant part of V can actually be 
ignored i.e. as far as the background expansion is con- 
cerned, a theory with a small scalar field value such as 
f{R) should be very close to simply having a cosmologi- 
cal constant [1^. A possible exception is if the potential 
violates dV/dip ~ SV/ip. 

In the non-relativistic (no time derivatives, small pe- 
culiar velocities) plus sub-Hubble hmit, the first order 
Einstein equations reduce to the analog of Eq. ([57)1 : 



V2($-v|) = 0. 



(E5) 



Writing ip = (p + Sip, the perturbed (though not neces- 
sarily first order) scalar field equation in the quasi-static 
limit is 



-irW^Stp = —- + aSnGSp . 
ooip 



(E6) 



analogous to Eq. 

To implement the surface integral approach to the 
problem of motion, we need some generalized concept 
of conservation in comoving space. Let us make the fol- 
lowing split: Gf/ = G^"" + 6Gt,\ where G^" is the FRW 
background piece. The perturbed (not necessarily small) 
piece 5Gp can be written as the sum of the first or- 
der part 5&^\ plus the rest. It can be shown that in 
the quasi-static, sub-Hubble limit, we have the following 
identities: 



- 
Go 



1 fc 



(E7) 



where the first expression assumes the quasi-static/non- 
relativistic, sub-Hubble limit while the second expression 
is exact. 

The perturbed Einstein equation, keeping all the non- 
linearities but subtracting out the background FRW 
piece, can be written as 



5t, 



ST, 



SttG 



(JG/ - <5G(i)/) . (E8) 



We therefore have the following modified 'conservation' 
laws of the pseudo-energy-momentum: 



- 
To 



1 ^ k 

-Tk 
3 



(E9) 



Note that indexes on Sifj, are raised and lowered by the 
background FRW metric. 

The mass of our object of interest is 



M 



(ElO) 



where we assume the mass is dominated by the overden- 
sity; the background contributes very little over our vol- 
ume. Note that x is comoving. The first conservation law 
can be used to show that M is nearly constant — as long 
as there is negligible fiux through the surface bounding 
our volume. Define then the center of mass coordinate 



d?xa^5tQx'/M , 



(Ell) 



X' = - / d^xa^Sto /M^ / d''xa^6ti/M, (E12) 



and the same conservation law tells us 
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where the ' denotes derivative with respect to conformal 
time. Finally, we can employ the second 'conservation' 
law to show 



a 

M 



-dj5t 



- 



- k " 

Tk 



(E13) 

The first term on the right can be rewritten as a surface 
integral using Gauss law. The second term, upon split- 
ting <I> = <I>o + $1 as usual, gives us something like adi^o 
times the mass from the background, which is negligible 
compared with the mass of the object. 

Carrying out the surface integral, we obtain the analog 



ofEq. 

M 



X* H X* 



^Af [1 + 2ea2](9j<lo , (E14) 



which holds in the quasi-static, sub-Hubble limit. The 
only difference from our earlier result is the appearance 
of a Hubble drag term, as expected. 

Transforming to Jordan frame, the equation is 



M 



^ — x"- 



-Ma,$o + (1 - e)aMa(po (E15) 



which is analogous to Eq. (|C16|) . 
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